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ABSTRACT 

This  study  was  conducted  at  the  Naval  Postgraduate  School 
to  investigate  the  formation  of  particle  flow  cells  at  mini- 
mum fluidization  gas  flow  rates  in  a  rectangular  gas-f luidized 
bed.   The  primary  objective  is  to  determine  if  multiple  cell 
formation  occurs  as  the  bed  width  to  height  ratio  increases . 
The  settled  bed  height  is  maintained  at  a  constant  level  for 
all  runs.   The  secondary  objective  is  to  determine  if  the 
presence  of  external  plate  heaters  at  the  end  walls  of  the 
particle  bed  will  alter  the  particle  flow  cell  formation 
patterns  when  energized  and  to  what  extent  they  affect  the 
fluid  mixing  (heat  transfer)  flow  patterns  within  the  par- 
ticle flow.   Pressure  measurements  across  the  bed  are  used 
to  determine  minimum  fluidization  points  and  thermocouple 
probe  measurements  determined  the  thermal  gradients  within 
the  bed  for  each  bed  configuration.   Visual  observations  on 
the  external  bed  movement  are  included  in  the  determination 
of  the  results  of  this  experimental  study. 
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NOMENCLATURE 

English  Letter  Symbols 

a  Constant,  dependent  on  tube  orientation  -  Reference  3  0 

A  Cross-sectional  area  of  bed  container, 
c 

A  Contact  area  of  packet  and  surface  -  Figure  18. 

B  Shape  factor  -  Reference  12. 

3i  Biot  number. 

C  Specific  heat  of  packet  or  quiescent  bed. 

C  f  Specific  heat  of  fluid. 

C  Specific  heat  of  solid  particle. 

CR  Constant,  correction  for  location  of  immersed  body  - 
Reference  34. 

C  Specific  heat  of  solid  particle. 

dp  Diameter  of  sphere  of  same  specific  area  as  that 
of  the  bed  particle. 

D  Diameter  of  the  fluidized  bed. 

D1  Effective  diameter  of  free  area  across  the  bed. 

DH  Diameter  of  an  immersed  heater. 

D.  Particle  diameter  of  sieved  fraction. 

D  Particle  diameter. 

D  Average  particle  diameter. 

D  Distance  between  two  successive  layers  -  Figure  16 . 

DT  Diameter  of  immersed  tube. 

D2Q  Diameter  of  immersed  tube  being  equal  to  20  mm. 

g  Acceleration  of  gravity. 
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G  Fluid  mass  velocity. 

G'  Fluid  mass  velocity  based  on  voidage  area. 

Ga  Galileo  number  (Archimedes  number) . 

G  Hypothetical  fluid  mass  velocity  required  to 
merely  expand  a  bed  of  particles . 

Gf  Fluid  mass  velocity  through  a  fluidizing  bed. 

G  -  Fluid  mass  velocity  for  minimum  f luidization . 

G  Optimum  fluid  mass  velocity  giving  a  maximum 

"  value  of  the  heat  transfer  coefficient. 

G  Fluid  mass  velocity  based  on  kinematic  viscosity. 

h  Heat  transfer  coefficient  for  the  contact 
resistance  at  the  wall. 

h  Empirical  constant  -  Reference  44. 

o  c 

h  Heat  transfer  coefficient  between  the  fluidized 
bed  and  the  surface. 

h  .  _  Instantaneous  local  heat  transfer  coefficient. 

Wlli 

h  Local  average  heat  transfer  coefficient. 

h  Maximum  heat  transfer  coefficient  between 
fluidized  bed  and  surface. 

H  Superficial  bubble-to-cloud  heat  transfer  coefficient 

k^  Fluid  thermal  conductivity. 

k  Thermal  conductivity  of  packet  or  quiescent  bed. 

k  Thermal  conductivity  of  particle. 

k  Thermal  conductivity  of  sphere. 

sp  J      c 

l  Distance  from  bed  top  in  packet  theory  -  Reference  47 

L  Fluidized  bed  height. 

L  Expanded  bed  height. 

LT  Length  of  heat  exchange  surface. 
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L  j.  Bed  height  at  point  of  minimum  f  luidization . 

m  Slope  -  Figure  12. 

m'  Slope  -  Figure  15. 

N  Fluidization  number. 

Nu  Nusselt  number. 

Nu  Maximum  value  of  Nusselt  number, 
max 

p  Empirical  constant  -  Reference  44. 

P  Center-to-center  distance  (pitch)  of  adjacent  tubes. 

Pr  Prandtl  number. 

q  Heat  flux. 

q  Heat  transferred  into  the  core  portion  of  the  bed  by 
an  interchange  of  particles  -  Figure  23. 

q  Heat  flow  rate  from  the  wall  into  a  packet. 

q  Instantaneous  heat  transfer  rate. 

q  Heat  taken  by  particles  flowing  parallel  to  the 

surface  in  a  second  zone  of  emulsion  of  thickness 
6   -  Figure  23. 

r  Fluctuation  ratio  -  highest  to  lowest  bed  level  in  a 
gas-f luidized  charge. 

r, Radius  of  sphere. 

sp  r 

R  Bed  expansion  ratio. 

Re  Reynolds  number. 

Ref  Reynolds  number  based  on  particle  free  fall  velocity 

Re  Modified  Reynolds  number  -  Reference  27. 

Re  -  Reynolds  number  at  minimum  fluidization. 

Re  Optimum  value  of  Reynolds  number  giving  a  maximum 

^  heat  transfer  coefficient. 

Re,  Reynolds  number  based  on  bed  diameter. 
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s  Average  replacement  of  packets  at  wall  by  side  mixing 

St  Stanton  number. 

t  Time. 

T,  Bulk  bed  constant  temperature. 

Tf  Fluid  film  temperature. 

T.  Gas  temperature  into  bed. 

T  Gas  temperature  out  of  bed. 

T  Temperature  of  wall  or  heater. 

w        r 

u.  Initial  fluid  velocity  with  no  particle  motion, 

u  Fluid  velocity  required  for  initial  expansion, 

u  r  Fluid  velocity  required  for  minimum  f luidization . 

u  Superficial  velocity  at  f luidization  conditions. 

u  Particle  velocity. 

P 

u  Fluid  velocity  during  f luidization. 

u  Particle  velocity  of  Packet  Model, 
s 

u  Terminal  fluid  velocity. 

V  Volume  displaced  by  particles. 

W  Bed  width. 

x  Axis  direction  -  Figures  21  and  22. 

x.  Weight  fraction  of  sieved  particles  for  D.. 

X  Packet  width  -  Figures  21  and  22 . 

y  Horizontal  pitch  of  tube  bundle  -  Reference  47. 

z  Vertical  pitch  of  tube  bundle  -  Reference  47. 


Greek  Letter  Symbols 

a      Constant  -  Reference  9 
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Y  Constant  -  Reference  9 . 

yf  Specific  weight  of  fluid. 

5  Boundary  layer  thickness  at  point  between  two  particle 
layers  -  Figure  16. 

r+  Average  boundary  layer  thickness  at  point  between  two 
points  -  Figure  16. 

6  Thickness  of  emulsion  layer. 

6_  Thickness  of  gas  layer. 

A  Minimum  distance  between  surfaces  of  adjacent  tubes 
(gap)  . 

Ap  Pressure  drop  across  bed  at  u-,. 

Ap  Pressure  drop  across  bed  at  u  . 

Ap  _  Pressure  drop  across  bed  at  u  - . 

Ap  Ideal  pressure  drop  across  the  bed  calculated  from  the 
buoyant  weight  of  the  bed. 

e  Void  fraction  at  f luidization  conditions . 

e.  Initial  particle  voidage  as  filled. 

e  Particle  voidage  during  expansion. 

e  -  Particle  voidage  at  minimum  f luidization . 

£Q  Void  fraction  at  f luidization  conditions. 

1  Fluidization  efficiency. 

9  Non-dimensional  absolute  bed  temperature. 

M  Fluid  viscosity. 

v  -  Fluid  kinematic  viscosity  at  minimum  fluidization. 

T"  Constant. 

P-r  Fluid  density. 

P  Density  of  packet. 

Pmf  Fluid  density  at  minimum  fluidization. 
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p  Particle  density. 

p  ,  Nonfluidized  bed  density. 

sb  z 

tj>  Sphericity  of  solid  particle 
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I.   INTRODUCTION 

If  one  considers  a  fluid  flowing  through  a  bed  of  solid 
particles,  several  forms  of  today's  fluidized  bed  technology 
are  created.   When  the  particles  are  stationary,  the  process 
is  called  flow  through  a  porous  medium  or  fixed  bed  fluidiza- 
tion.   If  the  fluid  velocity  is  such  that  the  particles 
become  entrained  by  the  fluid  and  are  removed  from  the  bed, 
the  process  is  then  called  pneumatic  conveying  or  entrained 
flow  f luidization.   At  some  point  between  these  two  conditions, 
the  particles  in  mixture  with  the  fluid  are  transformed  into 
a  fluid-like  state  causing  continuous  interaction  between  the 
two  mediums.   The  bed  particles  in  this  state  do  not,  however, 
leave  the  bed  container.   This  state  is  called  fluidized  bed 
f luidization. 

The  technique  of  fluidization  (all  forms)  is  used  for 
solid-fluid  interaction  processes  whenever  high  rates  of  heat 
and  mass  transfer  between  the  two  constituents  is  required  or 
high  rates  of  particle  transport  to  and  from  the  bed  is 
necessary.   The  advantages  in  the  use  of  fluidization  can  be 
listed  as  follows:  [6,12]  1 

(a)   The  formation  of  an  extremely  large  area  of  contact 
between  the  particles  and  fluid.   This  very  large  surface 


lumbers  in  brackets  indicate  references 
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area  permits  the  achievement  of  high  overall  rates  of  heat 
and  mass  transfer  between  the  particles  and  the  fluid. 

(b)  The  comparative  ease  with  which  the  fluidized 
particles  can  be  handled. 

(c)  The  reduction  of  temperature  gradients  to  negligible 
proportions  throughout  the  bulk  of  the  bed  as  a  result  of  the 
high  degree  of  particle  mixing  that  can  occur  in  gas-f luidized 
systems,  i.e.  a  very  high  effective  internal  thermal  conduc- 
tivity. 

(d)  The  occurrence  of  high  rates  of  heat  transfer  between 
the  fluidized  particles  and  an  immersed  surface. 

(e)  The  value  of  the  high  thermal  inertia  of  the  solids 
in  heat  transfer  operations  within  the  bed  as  a  whole. 

(f)  The  face  that  the  gas-f luidized  systems  approximate 
a  fluid  which  has  a  low  vapor  pressure  even  at  very  high  tem- 
peratures -  as  a  thermodynamic  fluid  it  lies  roughly  between 
gases  and  liquids. 

Inherent  with  the  fluidization  techniques  are  also  some 
limitations:   [6,12] 

(a)  The  process  demands  the  expenditure  of  power  for 
fluidization. 

(b)  Some  solids  cannot  be  fluidized  because  of  either 
an  excessive  tendency  for  particle  attrition  or  agglomera- 
tion to  occur.   The  first  will  produce  too  many  fines  by 
degradation  of  the  particles  which  will  be  lost  from  the  bed. 
The  latter  will  give  "agglomerated"  particles  too  large  to 
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fluidize.   Also,  there  are  those  materials  of  such  irregular 
shape  that  they  can  only  be  handled  by  spouting  bed  techniques. 

(c)  Operating  rates  are  limited  to  within  the  range  over 
which  the  bed  may  be  fluidized.   If  the  fluid  flow  is  too  low, 
the  bed  settles;  if,  on  the  other  hand,  the  rate  is  too  high, 
there  is  an  excessive  loss  of  particles  from  the  bed  by  carry- 
over . 

(d)  There  is  usually  a  limit  to  the  size  of  the  particles 
which  can  be  fluidized.   This  limit  is  approximately  6 . 5  mm 
(0.25  in)  down  to  a  few  microns  (0.001  in) .   It  is  particularly 
dependent  on  the  fluid  viscosity  and  particle/fluid  density 
ratio  as  well  as  the  fluid  flow  rate. 

(e)  True  counter  flow  operation  is  unobtainable  because 
of  the  high  degree  of  mixing  that  usually  occurs  within  a 
fluidized  bed. 

(f)  The  general  absence  of  thermal  gradients  throughout 
a  well  mixed  fluidized  bed  is  a  disadvantage  in  those  opera- 
tions where  thermal  gradients  are  required. 

(g)  The  dynamics  of  fluidized  systems  are  not  sufficiently 
well  understood  to  permit  confident  scale-up  to  large  scale 
units  from  small  scale  test  data  when  reaction  conditions  are 
critical. 

Up  until  this  time  the  term  fluid  has  been  used  in  describ- 
ing the  transport  constituent  in  f luidization.   To  be  more 
specific  there  are  two  fluids  which  have  received  widespread 
use  in  fluidized  bed  technology,  air  (gas)  and  water  (liquid) . 
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The  basic  differences  are  in  the  dynamics  and  the  mass  and 
heat  transfer  of  the  solid  particles  suspended  in  the  gas- 
fluidized  or  liquid-f luidized  bed.   This  is  primarily  a  result 
of  the  difference  in  the  particle  to  fluid  density  ratios.   In 
the  gas-f luidized  bed  the  ratio  will  usually  be  on  the  order  of 
a  thousand  while  in  a  liquid-f luidized  bed  the  particles  are 
only  several  times  heavier  than  the  fluid. 

Additionally,  the  fluid  flow  characteristics  will  affect 
the  overall  mass  and  heat  transfer  bed  characteristics.   In 
the  liquid-f luidized  system,  the  fluidization  velocities  are 
usually  of  such  a  magnitude  that  the  fluid  inertial  terms  are 
considered  insignificant.   This  is  not  the  case  with  a  gas- 
fluidized  system,  especially  when  large  particle  sizes  are 
used.   As  the  velocities  increase  the  inertial  fluid  losses 
and  form  drag  become  important  in  determining  the  bed  charac- 
teristics . 

With  the  smooth  manner  in  which  liquid-f luidized  beds 
operate,  i.e.  no  gross  instabilities  in  flow,  the  bed  is 
sometimes  referred  to  as  a  homogeneously  fluidized  bed  or 
a  particulately  fluidized  bed.   The  gas-f luidized  beds  on 
the  other  hand  usually  have  unstable  flow  above  some  minimum 
fluidization  condition,  characterized  by  violent  agitation 
of  the  particles  accompanied  by  bubbling  and  channeling  of 
the  gas.   As  a  result,  the  gas-f luidized  bed  is  sometimes 
called  an  aggregate  or  bubbling  fluidized  bed. 
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The  bulk  of  recent  research  has  been  in  the  area  of 
gas-f luidized  beds  vice  liquid-fluidized  beds.   Whether  this 
trend  will  shift  remains  to  be  seen.   Nevertheless,  this 
investigation  will  be  in  the  gas-f luidized  field  and  all 
further  discussions  will  be  based  on  available  gas-solid 
research  data. 

In  a  gas-f luidized  bed  the  particles  are  supported  above 
a  gas  distributor  plate  which  permits  the  gas  to  enter  and 
flow  upward  through  the  bed.   As  the  gas  flow  increases, 
there  will  be  a  pressure  drop  across  the  bed  which  will 
eventually  reach  a  point  when  it  just  balances  the  downward 
gravitational  force  of  the  particle  bed.   When  this  point  is 
reached,  the  bed  will  expand  and  the  particles  separate. 
This  is  referred  to  as  the  minimum  fluidization  point.   Fur- 
ther gas  flow  rate  increases  will  result  in  bubble  formation 
which  will  rise  through  the  bed  and  burst  at  the  surface  not 
unlike  a  vapor  bubble  in  a  boiling  liquid. 

The  separation  of  the  particles  creates  a  large  particle 
surface  in  direct  contact  with  the  gas  and  the  bubble  flow 
induces  vigorous  and  thorough  mixing  of  the  particles.   This 
high  gas-particle  contact  area  results  in  the  high  heat 
transfer  rates  which  have  been  experienced.   Because  the 
particle  flow  takes  place  in  the  fluidizing  gas,  a  secondary 
fluid-mixing  phenomenon  will  be  induced  through  the  particle 
flow  pattern.   The  extent  of  this  secondary  mixing  and  its 
effect  on  heat  transfer  will  depend  on  whether  the  particles 
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are  dispersed  in  the  bed  more  or  less  as  individual  particles 
or  clumps  of  particles.   As  an  individual  particle,  disper- 
sion by  the  fluid  mixing  would  be  relatively  small.   As  an 
agglomerate  of  particles,  however,  the  result  is  a  dragging 
of  gas  pockets  through  the  bed  resulting  in  perturbations  to 
the  particle  flow  patterns.   Thus,  the  extent  of  fluid  mixing 
is  related  to  the  mode  and  quality  of  f luidization.   This  is 
an  area  that  has  been  scarcely  investigated.   Two  properties, 
though,  that  have  been  established  as  a  result  of  fluid 
mixing  investigations  are  that: 

(a)  Owing  to  fluid  back-mixing,  temperature  driving 
forces  are  reduced. 

(b)  Analogously  to  heat  transfer,  driving  forces  are 
also  reduced  with  mass  transfer. 

Based  on  previous  research  the  normal  particle  flow 
pattern  in  small  diameter  cylindrical  gas-f luidized  beds  is 
as  shown  in  Figure  1.   For  gas  flow  rates  which  are  signif- 
icantly greater  than  the  minimum  fluidization  flow  rate  the 
solids  tend  to  move  downward  at  the  bed  walls  and  the  passage 
of  the  individual  particle  is  readily  observed.   The  particles 
will  move  down  the  wall  and  then  submerge  into  the  body  of 
the  bed.   Around  this  entry  point  other  particles  will  emerge 
from  the  bed  and  travel  down  the  wall  until  they  also  merge 
back  into  the  bed.   This  process  continues  until  the  bottom 
of  the  bed  is  approached.   At  some  point  near  the  bottom  of 
the  bed  an  inactive  or  dead  region  will  form  in  the  particles. 
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The  size  of  this  region  will  depend  on  the  height  of  the  bed, 
fluid  flow  rate  and  bed  container  geometry.   In  deep  beds 
the  flow  regions  are  well  defined  while  in  very  shallow  beds 
the  downward  flow  pattern  is  much  less  clearly  defined  with 
the  bed  appearing  to  be  more  turbulent  and  in  a  homogenous 
state  of  agitation. 

In  the  literature  reviewed,  when  the  term  shallow  bed  was 
used  it  generally  was  referring  to  a  particle  bed  whose  depth 
was  on  the  order  of  one  and  a  half  feet  or  less  depending  on 
the  fluid  inlet  cross-sectional  area.   In  general,  the  major 
features  distinguishing  shallow  beds  from  deep  beds  can  be 
described  as:   [55] 

(a)  Violent  coalescence  of  bubbles  taking  place  through 
the  majority  of  the  shallow  bed  depth. 

(b)  The  jetting  region  next  to  the  distributor  plate 
occupying  a  significant  proportion  of  the  shallow  bed  depth. 

(c)  The  thermal  "entry  region"  occupying  a  significant 
proportion  of  the  shallow  bed  depth. 

These  features  would  suggest  that  the  particle  motion 
within  the  shallow  bed  may  well  be  significantly  greater 
than  in  a  deep  bed  and  that  the  bubble  growth  will  largely 
be  completed  within  a  very  small  proportion  of  the  bed  depth. 
Thus  the  residence  times  of  the  particles  at  the  heat  trans- 
fer surface  would  be  considerably  shorter  than  in  a  deep  bed. 

In  the  flow  description  of  deep  beds,  the  bed  particles 
will  move  upward  in  the  center  of  the  bed  either  in  a  spouting 
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type  pattern  or  in  a  multi-channel  pattern.   It  should  be 
noted  that  if  the  bed  vessel  is  operated  in  a  configuration 
such  that  its  containment  walls  are  not  parallel  to  the 
fluid  distributor  plate  inlet  gas  flow  direction,  then  the 
particle  flow  pattern  will  be  adversely  affected.   Only  a 
few  degrees  of  deviation  from  this  parallel  configuration 
will  lead  to  particle  flow  pattern  short-circuiting.   If  the 
operating  gas  velocities  are  further  increased  while  this 
configuration  is  present,  then  slugging  will  initialize. 

As  was  previously  stated,  fluid  mixing  in  gas-f luidized 
beds  will  occur  as  a  result  of  the  particle  mixing  at  the 
onset  of  f luidization.   Unfortunately,  the  gas  does  not 
circulate  in  the  smooth  manner  shown  in  Figure  2.   Instead, 
there  will  be  many  eddies  within  the  central  core  of  the 
particle  flow  pattern  as  well  as  horizontal  components 
radiating  out  from  the  core.   Additionally,  part  of  the  gas 
may  rise  rapidly  in  the  form  of  bubbles,  rising  ahead  of  the 
slower  moving  interstitial  gas  components.   Therefore, 
Figure  2  should  only  be  viewed  as  an  ideal  flow  pattern  at 
minimum  f luidization  for  the  bed. 

To  date,  particle  flow  patterns  with  the  resulting  par- 
ticle and  fluid  mixing  criteria  have  been  explored  only  in 
equipment  of  relatively  small  diameter  (less  than  one  foot) . 
The  experimental  apparatuses  that  were  used  were  circular 
in  nature  with  either  no  external  heat  source  or  a  submerged 
cylindrical  (tube)  heat  source.   The  data  obtained  was 
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Figure  1.   Observed  particle  circulation  pattern 


Gas    Out 

t    t  t  t    1 


\f/ 


+     i 


l     * 


//n_\\\ 


Gas    In 


Fluid 

.Ci  rcu  1  at  i  on 
Pattern 


.Di  str  ibutor 
Plate 


Figure  2.   Greatly  simplified  gas  flow  pattern 
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usually  a  byproduct  of  the  investigation  to  determine  heat 
transfer  coefficients  for  fluidized  beds  with  various  heat 
sources  and  bed  geometries . 

This  thesis  investigation  was  conducted  to  determine 
the  particle  flow  patterns  developed  in  a  gas-f luidized  bed 
of  variable  rectangular  geometry.   Secondly,  to  determine 
what  effect  an  external  plate  heat  source  has  on  the  particle 
and  fluid  flow  patterns  developed  within  the  rectangular  bed. 
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II.   OBJECTIVES  AND  METHOD  OF  APPROACH 

A.  GENERAL  STATEMENT  OF  THE  PROBLEM 

The  specific  problem  may  be  defined  as  an  investigation 
into  the  formation  of  particle  flow  cells  at  minimum  gas 
fluidization  flow  rates  in  a  rectangular  gas-f luidized  bed. 
The  primary  objective  is  to  determine  if  multiple  cell  for- 
mation occurs  as  the  bed  width  to  height  (W/L  )  ratio 
increases.   The  secondary  objective  is  to  determine  if  the 
presence  of  plate  heaters  at  the  end  walls  of  the  particle 
bed  will  alter  the  particle  flow  cell  formation  patterns 
when  energized  and  to  what  extent  they  affect  the  fluid 
mixing  (heat  transfer)  flow  patterns  within  the  particle 
flow. 

B.  METHOD  OF  APPROACH 

The  methods  used  for  this  investigation  were  experimental 
and  analytical.   Using  available  research  literature  a  deter- 
mination was  made  of  the  minimum  fluidizing  velocity  (u  f) 
for  the  particular  f luidized  bed  configuration.   This  value 
of  u -  was  based  on  recorded  and  analytical  pressure  drops 
across  the  distributor  plate  and  bed  in  addition  to  the 
expanded  bed  heights.   Visual  confirmation  of  the  minimum 
fluidizing  conditions  was  then  compared  with  the  analytical 
estimates  of  the  required  pressure  drops,  bed  height  and  air 
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velocities.   Once  an  initial  particle  flow  cell  was  estab- 
lished, the  bed  geometry  and/or  activation  of  the  plate 
heaters  would  be  used  in  an  attempt  to  alter  its  configura- 
tion.  A  grid  of  variable  depth  copper-constantan  thermo- 
couples was  used  to  map  the  thermal  distribution  within  the 
particle  flow  cell  to  determine  the  fluid  mixing  properties 
of  the  cells. 

This  thesis  will  not  be  involved  with  the  heat  transfer 
coefficient  determination  so  prevalent  in  recent  experimenta- 
tion; rather  it  will  focus  only  on  how  the  particle  flow  cells 
relate  to  the  bed  geometry  chosen  and  how  the  resulting  fluid 
mixing  affects  the  heat  transfer  within  the  bed.   It  should 
be  noted  that  direct  measurement  of  the  solid  particle  temper- 
atures is  very  unreliable  and  that  the  temperatures  recorded 
will  be  considered  to  be  that  of  the  air  at  a  particular 
point  within  the  grid  of  the  bed.   Because  of  the  difference 
in  the  particle  to  fluid  density,  however,  the  two  will  be 
considered  to  be  at  equilibrium  at  the  point  measured  for 
this  experiment. 
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III.   SUMMARY  OF  PREVIOUS  RESEARCH  IN 
GAS-FLUIDIZED  BED  TECHNOLOGY 

A.   GENERAL  DESCRIPTION  OF  GAS-FLUIDIZED  BED  BEHAVIOR 
1 .   The  Phenomenon  of  Fluidization 

When  the  bed  configuration  is  initially  filled  with 
particles,  an  initial  voidage  (e^)  will  result  between  the 
particles.   As  the  fluidizing  fluid,  a  gas  in  this  case,  is 
initiated  at  some  superficial  velocity  (u,)  there  will  be  no 
movement  of  the  bed.   However,  a  pressure  drop  (Ap,)  will 
occur  and  gradually  increase  at  a  linear  rate.   The  magnitude 
of  this  pressure  drop  will  be  determined  by  the  gas  flow  rate 
and  the  characteristics  of  the  bed. 

Eventually  a  condition  will  be  reached  when  the  gas 
flow  rate  (u  )  has  increased  to  a  point  at  which  expansion 
of  the  bed  initiates.   The  voidage  will  increase  to  some  e 
and  the  pressure  drop  can  be  calculated  from  the  equation 

Ape  =  Le  (  1  -  ee  )  (  Ps  "  Pf  )  (3.1) 

(Note  -  The  units  to  be  used  with  all  text  equations  can  be 

found  in  Appendix  D.)   While  the  bed  is  expanding,  there 

still  will  have  not  been  any  particle  motion  taking  place. 

As  the  gas  flow  rate  is  increased  further,  a  point 

will  next  be  reached  when  the  voidage  becomes  e  ,-  and  the 
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particles  are  at  the  point  of  fluidization  onset.   The 
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velocity  is  called  the  minimum  fluidization  velocity  (u  f) 
and  is  critical  to  most  bed  behavior  calculations.   Further 
gas  velocity  increases  (uQ)  will  fluidize  the  bed  causing 
further  voidage  increases  (e„)  until  a  terminal  velocity 
(u. )  is  reached  as  defined  by  the  bed  container  geometry  or 
particle  entrainment  takes  place.   By  the  time  u.  has  been 
reached  the  particle  motion  will  be  intense  although  the 
pressure  drop  across  the  bed  will  remain  at  the  level  of 
minimum  fluidization  inception.   This  condition  is  referred 
to  as  dense-phase  turbulent  fluidization  at  its  inception 
(un)  but  degenerates  to  a  dilute-phase  when  the  gas  rate  is 
such  that  the  particles  become  entrained  in  the  gas.   On  the 
average  the  ratio  of  the  terminal  velocity  (u, )  to  the 
minimum  fluidization  velocity  (u ~)  is  about  seventy  (70) . 
The  terminal  velocity  is  also  a  function  of  the  particle 
diameter,  shape  and  density  as  well  as  the  properties  of  the 
fluidizing  fluid. 

Parent  et  al.  [52]  has  characterized  three  types  of 
behavior  with  the  gas-solid  mixtures  in  the  dense-phase 
fluidization  region.   Under  certain  conditions,  usually  for 
small  particle  sizes,  the  gas  tends  to  channel  through  the 
bed  without  causing  movement  of  the  bed.   Under  other  condi- 
tions, such  as  large  bed  heights,  large  particle  sizes  and 
high  gas  flow  rates,  slugging  takes  place  which  is  analogous 
to  the  turbulence  observed  in  a  boiling  liquid.   The  third 
type  of  behavior  is  described  as  the  mechanically  smooth 
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f luidization.   In  this  case  the  particles  move  along  twisting, 
ever-changing  paths.   Occasionallly  gas  bubbles  appear  in  the 
bed  but  there  is  no  rough  bumping.   It  is  this  case  that  the 
pressure  drop  across  the  bed  is  nearly  equal  to  the  suspension 
head  opposing  the  flow  (drag) .   A  word  of  caution  is  warranted, 
however,  in  that  in  research  conducted  by  Lewis  et  al .  [42] 
the  conclusion  that  the  pressure  drop  is  essentially  equal  to 
the  weight  of  the  particle  bed  would  appear  to  be  satisfactory 
only  for  units  with  small  L/D  ratios.   It  appears  to  be  in 
considerable  error  for  large  L/D  values  or  for  beds  having 
large  static  effects  present. 

Here  again  several  researchers  differ  in  their  opin- 
ion as  to  the  cause  of  this  phenomenon.   While  Lewis  et  al . 
[42]  explained  it  as  an  indication  of  frictional  drag  of  the 
solid  particles  on  the  walls  of  the  unit,  Toomey  and  Johnstone 
[62]  attribute  it  to  kinetic  energy  losses  caused  by  the  colli- 
sion between  the  particles  and  Wilhelm  and  Kwauk  [76]  con- 
cluded it  was  caused  by  forces  which  are  electrostatic  and/or 
fluid  dynamic  in  origin.   Since  this  thesis  is  operating  in 
the  shallow  bed  region  it  is  not  anticipated  that  this  excess 
pressure  drop  will  be  of  concern. 

2 .   Voidage  at  Minimum  Fluidization 

The  first  bed  characteristics  that  must  be  determined 
for  an  understanding  of  gas-f luidized  bed  operation  is  that 
of  the  bed  voidage  at  minimum  fluidization  (e  ^) .   By  defini- 
tion e  -  is  that  particle  separation  to  which  the  bed  must 
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expand  before  particle  motion,  as  induced  by  fluid  flow,  can 
initialize.   There  are  two  methods  which  can  be  used  to  deter- 
mine this  factor. 

The  first  method  is  based  on  a  collection  of  data  from 
several  researchers  that  has  been  graphically  combined  by  Leva 
[36]  and  reproduced  here  in  Figure  3.   For  a  given  particle 
diameter  (D  )  and  composition  a  value  for  the  minimum  fluidiza- 

tion  voidage  (e_^)  can  be  obtained.   The  obvious  drawback  for 
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this  method  is  that  it  depends  solely  on  the  material  chosen 
for  the  particles  being  one  that  has  had  previous  experimenta- 
tion accomplished. 

The  second  method  is  used  with  those  beds  for  which 
particle  z    -  data  is  not  available.   It  requires  that  the  bed 
be  subjected  to  the  fluidizing  gas  until  particle  motion 

initializes,  at  which  time  the  bed  height  (L,-)  is  recorded. 

3     mr 

The  minimum  fluidization  voidage  (e  ^)  can  then  be  calculated 

J     mr 

from  the  equation 

v   Y* 
Emf  -  1 •  (3.2) 

Lmf  Ac  (  ps  "  pf  } 

The  accuracy  of  this  estimate  will  depend  on  the  mass  charac- 
teristics of  the  bed  particles,  i.e.  sphericity,  whether  solid 
or  hollow,  degree  of  particle  porosity,  etc.   Particles  of 
high  porosity  are  usually  unsuitable  for  fluidized  bed  use 
as  they  have  a  tendency  to  break  up  with  increased  fluid  flow. 
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Figure  3.   Values  o-f  £mf  vs .  particle  diameter. 


34 


Once  a  value  for  the  minimum  fluidization  voidage 
has  been  estimated,  additional  bed  behavior  characteristics 
can  be  predicted  and  the  bed's  suitability  for  the  intended 
process  determined. 

3 .   Pressure  Drop 

Another  excellent  indication  of  fluidization  behavior 
of  the  particle  material  may  be  obtained  from  the  pressure 
drop-flow  relationship  of  the  expanded  bed  [36,40,64].   Ideally, 
the  pressure  drop  across  the  bed  should  increase  linearly  until 

the  voidaqe  has  reached  a  maximum  value  (e    )  at  which  time  the 
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pressure  differential  would  level  off  for  the  remaining  phases 
of  fluidization.   See  Figure  4. 

Unfortunately,  the  idealization  never  occurs  as  a 
result  of  three  major  factors;  namely  channeling,  slugging 
and  distributor  design.   With  channeling  a  pressure  drop 
occurs  (Figure  5)  when  £   is  reached.   The  magnitude  of  this 
drop  will  depend  on  the  extent  of  the  channeling  phenomenon 
but  eventually  a  minimum  drop  point  is  reached  and  the  pres- 
sure differential  then  increases  to  some  new  lower  steady 
state  value . 

The  channeling  can  be  either  "through  channeling"  in 
which  it  extends  completely  through  the  bed  or  "intermittent 
channeling"  where  it  extends  through  part  of  the  bed.   For- 
tunately, as  the  gas  flow  rate  is  increased  the  tendency  for 
the  bed  to  have  channeling  decreases.   Some  factors  [38] 
which   tend  to  promote  channeling  are  poor  distributor  design, 
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Figure    4.       Ideal    pressure-drop-f  low    diagram. 
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Figure  5.   Pressure-drop-f  1  ou»  diagram  -For  channeling. 
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the  shape  of  the  bed  particles  (less  sphericity  the  greater 
the  channeling) ,  moisture  in  the  bed  and  the  rate  of  fluid 
flow  through  beds  of  small  particles.   If  the  bed  is  moist, 
the  particles  will  tend  to  clump  together,  and  the  gas  flow 
will  channel  between  the  clumps.   At  higher  flow  rates  the 
increased  agitation  of  the  bed  apparently  destroys  the 
channels  as  soon  as  they  are  formed. 

With  slugging  a  condition  exists  in  which  the  bubbles 
of  gas  combine  to  a  size  approaching  the  order  of  magnitude 
of  the  diameter  of  the  confining  vessel.   As  a  result,  the 
particle  layers  or  slugs  of  granular  particles  will  move 
upward  through  the  bed  like  a  piston  until  they  disintegrate 
at  the  bed  top.   When  this  happens,  the  pressure  drop  will 
increase  above  the  predicted  steady  state  value  (Ap  ) ,  see 
Figure  6.   If  the  flow  rate  is  increased  further,  the  pres- 
sure drop  rise  continues  and  the  bed  will  slug  even  more 
severely. 

The  third  factor,  namely  distributor  design,  was 
well  investigated  by  Grohse  [27]  who  stated  that  for  the  mul- 
tiorifice  plate  the  point  of  incipient  fluidization  is  not 
well  defined.   The  gradual  pressure-drop  increase  up  to  the 
theoretically  demanded  suggests  that  the  bottom  portion  of 
the  bed  is  not  fluidized  to  the  same  extent  as  the  upper 
region.   With  the  screen  distributor,  a  certain  amount  of 
channeling  probably  occurred.   There  is  the  possibility  that 
the  holes  in  the  screen  were  locally  blocked;  this  could 
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create  a  condition  favorable  for  channeling.   With  the  porous 
plate  the  point  of  incipient  fluidization  is  well  defined 
and  there  is  no  tendency  for  channeling.   See  Figure  7. 

Since  the  point  of  minimum  fluidization  is  the  key 
starting  point  on  fluidization  studies,  it  seems  extremely 
important  to  be  able  to  predict  that  point  both  theoretically 
and  through  experimental  results.   Therefore,  the  distributor 
design's  effect  on  the  pressure  drop  at  u  ^  becomes  a  very 
important  factor.   According  to  Kunii  and  Levenspiel  [34], 
distributors  should  have  a  sufficient  pressure  drop  to 
achieve  equal  flow  through  all  openings.   As  a  result,  the 
pressure  drop  across  the  distributor  plate  by  direct  measure- 
ment should  be  at  least  10%  of  the  pressure  drop  of  the  bed 
to  ensure  effective  fluidization.   Theoretically  the  pressure 
drop  at  the  minimum  fluidization  velocity  (Ap  J  can  be 

found  by  substituting  L  -  and  £  -  for  the  values  of  L   and 

■•  mr      mf  e 

e   in  equation  (3.1) . 

4 .   Minimum  Fluidization  Velocity 

The  ability  to  predict  the  point  of  initial  fluidi- 
zation reliably  is  of  basic  importance  in  virtually  all 
fluidized-process  studies  and  designs.   By  operating  at  mass 
velocities  which  are  slightly  higher  than  the  minimum  fluid- 
ization mass  velocity  (G  ^)  attrition  and  elutriation 
(stratification)  losses  will  be  minimized.   The  by-product 
of  the  low  flow  rates  is  then  a  low  particle  velocity  in  the 
bed  and  correspondingly  a  relatively  low  bed  wall  heat  transfer 
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The  general  equation  [36]  for  the  minimum  fluidization 
mass  flow  rate  is 

688  Dp1-82  {  pf  (  ps  -  Pf  )  }°-94 
G    = ,        (3.3) 

mr  0.8  8 


with  D  ,  the  average  particle  size, 


D   =  .  (3.4) 

p 

c  x. 

D. 

l 


Equation  (3.3)  can  only  be  used,  however,  when  the  Reynolds 
number 

0.03445  D   G  c 

p   mf 

Remf  =  — (3.5) 


is  less  than  ten.   Should  the  value  of  Re  f  be  greater  than 
ten,  then  the  alternative  is  to  use  the  nomograph  [72]  of 
Figure  8,  with  the  calculated  value  of  D  ,  y  and  Pf(ps  ~  Pf) 
to  recalculate  the  value  of  G  f.   If,  however,  the  value  of 
D,u   or  pf(p   ~  Pf)  does  not  appear  on  the  nomograph,  then 

IT 

go  to  Figure  9  [72]  and  with  the  calculated  value  of  Re  ,., 
from  equation  (3.5),  determine  the  correction  factor  to  be 
applied  to  the  value  of  G  ,  calculated  from  equation  (3.3) . 
The  new  value  of  G  f  is  then  calculated  from  the  equation 

mf   new      mf   correction     mf   calculated*    *  *  ' 
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NOTE-Left  (or  right)  scale  of  (ps-p.p)p.p 
corresponds  to  left  (or  right) 
scale  of  Gm.f 
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It  should  be  mentioned  that  there  are  numerous  corre- 
lations proposed  for  the  minimum  fluidization  mass  velocity 
(G  -)  that  are  applicable  to  the  experimental  results  of 
various  researchers.   These  predictions  are  not  in  very  close 
agreement  and  are  a  function  of  different  apparatus  construc- 
tions, fluids  used,  particles  used,  etc. 

One  variation  considered  worthy  of  note  is  a  result 
of  the  work  accomplished  by  Wen  and  Yu  [73].   They  demonstrated 
that  when  particles  of  more  than  two  different  sizes  are  fluid- 
ized,  they  either  separate  completely,  forming  distinct  layers 
or  intermix  forming  one  layer  of  rather  uniformly  distributed 
particles.   For  a  bed  containing  two  different  particle  sizes, 
if  the  particle  diameter  ratio  is  less  than  1.3,  the  particles 
were  observed  to  intermix.   If,  however,  the  particle  size 
ratio  was  greater  than  1.3,  the  smaller  particles  fluidized 
separately  above  the  larger  particles.   When  this  occurs, 
the  voidage  must  be  obtained  for  each  layer  and  special  pro- 
cedures, delineated  in  Reference  [73],  used  to  compute  the 
fluid  mass  velocity. 
5  .   Bed  Expansion 

This  area  of  bed  behavior  is  important  in  that  it 
resolves  the  minimum  bed  chamber  height  that  must  be  designed 
to  enable  fluidization  to  occur.   Unfortunately,  although 
gas-f luidization  processes  are  more  prevalent  than  liquid- 
fluidization  processes,  the  mechanics  of  the  gas-f luidization 
process  is  not  as  well  understood  in  comparison  to  the  liquid 
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process.   It  has  been  found  that  the  liquid  systems  will  ex- 
pand smoothly  almost  to  a  voidage  of  1.0.   With  a  gas  system, 
however,  somewhere  along  the  course  of  the  expansion,  a  point 
of  mechanical  instability  is  reached  which  degrades  the  expan- 
sion process.   The  location  of  this  trip  point  or  the  reason- 
ing as  to  why  it  occurs  has  so  far  eluded  researchers .   There 
is  some  speculation,  however,  that  it  may  be  caused  by  the 
formation  of  aggregates  of  particles  in  the  bed  and  that  it 
is  somehow  related  to  the  fluid  velocity. 

Based  on  analytic  and  experimental  data,  Leva  et  al . 
[36,40]  formulated  a  simple  relationship  between  the  fluid 
mass  velocity  and  bed  particle   voidage  that  enabled  the 
generation  of  a  prediction  model  for  bed  expansion.   Some 
examples  of  this  gas-solid  expansion  data  can  be  found  in 
Figure  10 .   What  is  shown  is  that  at  higher  flow  regions  the 
data  will  align  along  curves  rather  than  linear  lines  and  at 
intermediate  flow  regions  the  reverse  is  true.   Note  that  as 
the  voidage  increases,   (1  -  e)/e3  will  decrease.   Also,  at 
elevated  flow  rates,  the  slopes  (m)  become  less  negative  in- 
dicating that  aggregation  becomes  more  prevalent  while  the 
linear  lines  indicate  little  or  no  aggregation  (intermediate 
flow  region) .   It  is  also  noted  that  small-particle  beds  may 
be  expanded  over  substantially  wider  flow  ranges  than  larger- 
particle  beds  before  aggregation  conditions  will  set  in. 
This  is  indicated  by  the  fact  that,  with  smaller  particles, 

departure  from  the  straight  line  occurs  at  higher  values  of 
the  reduced  mass  velocity  than  with  larger  particles . 
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Figure    10.       Gas-particle    expansion    data 
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As  far  as  practical  methods  of  correlation  are  con- 
cerned, the  data  is  most  cases  may  be  represented  by  straight 
lines.   This  will  produce  the  relationship  between  slopes  (m) 
and  particle  size  (D  )  indicated  in  Figure  11  [36,40].   This 

cr 

correlation  is  based  on  virtually  all  of  the  available  bed- 
expansion  literature  data  to  date. 

To  calculate  a  predicted  expanded-bed  height  of  the 
gas-f luidized  particles,  two  basic  methods  may  be  used.   The 
choice  of  method  will  depend  on  the  bed-particle  diameter. 
However,  a  maximum  mean  particle  diameter  (D  )  of  about  0.015 
in  is  recommended,  otherwise  the  fluidization  efficiency  will 
be  so  low  as  to  promote  slugging.   Fluidization  efficiency 
(n)  is  defined  by  the  equation 

Gx  -  G 

f    e 

n  =  ,  (3.7) 

Gf 

where  Gf  and  G   are  respectively  defined  as  the  fluidizing 
gas  mass  velocity  flowing  through  a  bed  at  a  given  instant  in 
time  and  the  gas  mass  velocity  required  to  initially  expand 
the  bed.   Once  slugging  sets  in  the  bed  expansion  height 
will  begin  to  fluctuate  as  the  gas  bubbles  expand  and  collapse 
When  this  observed  condition  occurs  no  correlation  to  date 
has  been  found  to  give  a  predictable  bed  height, 
a.   Method  One 

(1)  With  D  known,  m  may  be  found  from  Figure  11. 

(2)  Determine  e  f  from  Figure  3  or  equation  (3.2) 
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(3)  Calculate  G  -  using  equations  (3.3)  and/or 
(3.6)  or  Figure  8. 

(4)  Draw  fluidization  line  on  log-log  paper  with 
slope  m,  see  Figure  12  [40]. 

(5)  For  the  chosen  value  of  Gf,  find  the  corre- 
sponding value  of  ef. 

(6)  Solve  the  voidage  equation,  equation  (3.2)  , 

for  L,.  by  substituting  L£  and  e-  for  L  ,  and  e  -. 
f  2  *      f      f      mf      mf 

b.   Method  Two 

(1)  Calculate  G  f  using  equations  (3.3)  and/or 
(3.6)  or  Figure  8. 

(2)  Calculate  the  reduced  mass  velocity,  Gf/G  -. 

(3)  Go  to  Figure  13  [36]  with  the  reduced  mass 

velocity  and  D   to  find  the  fluidization  efficiency  (n) . 

P 

(4)  Go  to  Figure  14  [36]  with  p  and  D   to  evaluate 
the  bed  expansion  ratio,  R. 

(5)  Calculate  the  static  bed  height  L  -  from 
equation  (3.2)  . 

(6)  Multiply  L  -  by  R  to  get  the  expanded  bed 
height,  Lf. 

Once  the  expanded  bed  height  has  been  calculated 
there  still  remains  one  check  to  be  performed  on  its  validity. 
This  is  a  check  on  the  extent  of  height  fluctuations  that  may 
occur  for  the  flow  rate  chosen.   It  has  been  found  that  under 
conditions  of  high  gas  flow  rates  that  an  oscillation  of  the 
expanded  bed  height  can  occur.   This  concept  has  been  termed 
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the  fluctuation  ratio,  r,  the  ratio  of  the  highest  to  lowest 
levels  which  the  bed  top  occupies  for  a  particular  gas  flow 
rate.   The  value  of  r  can  be  determined  from  the  equation 

m'  (  G.  -  G  -  ) 

r     mr 

r  =  exp   (3.8) 

mf 

where  slope  m'  is  related  to  the  particle  diameter  as  shown 
in  Figure  15  and  is  based  on  data  of  Lewis  et  al .  [42].   The 
closer  the  fluctuation  ratio  to  1.0  the  more  accurate  the 
expanded  bed  height  estimate. 

B.   HEAT  TRANSFER  PHENOMENON  MODELS 

There  is  little  agreement  in  the  reported  experimental 
values  of  the  heat  transfer  coefficients  in  the  literature. 
The  reasons  can  best  be  summarized  as  follows:   1)  a  variety 
of  contacting  patterns  between  particles  and  the  fluid  re- 
sulted in  different  models  from  which  to  base  the  average 
heat  transfer  coefficient;  2)  the  temperature  measurements 
taken  in  the  beds  by  thermocouples  were  interpreted  differ- 
ently (usually   the  readings  were  related  to  the  gas  temper- 
ature in  some  manner) ;  or  3)  that  direct  measurement  of 
solid  particle  temperatures  are  unreliable  with  present 
state-of-the-art  methods,  especially  for  the  smaller  size 
particles  (less  than  100  microns) . 

As  a  result,  due  to  the  generally  low  accuracy  of  in-bed 
temperature  measurements  and  the  oversimplifications  in  the 
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flow  models,  heat  transfer  coefficients  so  far  reported  have 
varied  widely.   Although  there  is  no  complete  theory  that 
can  predict  the  heat  transfer  data  for  fluidized  beds,  there 
are  several  theories  that  correlate  the  published  data.   Each 
resulting  heat  transfer  model  has  some  experimental  or  design 
situation  for  which  it  is  appropriate.   What  follows  is  an 
overview  of  these  various  models  all  of  which  are  based  on 
cylindrical  bed  geometries . 

1 .   Steady  State  Conduction  Across  the  Gas  Film 

In  this  model  heat  is  conducted  through  the  gas 
boundary  layer  near  the  heat  exchange  surface.   In  order  to 
explain  the  high  values  of  the  heat  transfer  rates  obtained 
experimentally,  the  gas  layer  is  assumed  to  be  scoured  by 
solids  moving  along  the  heat  exchange  surface,  thus  decreas- 
ing the  boundary  layer  thickness  near  the  wall.   This  model 
was  developed  by  Leva  et  al .  [3  7],  Dow  and  Jakob  [13],  and 
Levenspiel  and  Walton  [41] . 

Dow  and  Jakob  [13]  studied  the  gas-f luidization  mix- 
tures and  classified  them  as  either  a  dense  or  lean  phase 
mixture  based  on  the  gas  velocity.   The  dense  phase  occurred 
at  gas  velocities  close  to  u  ^  while  the  lean  phase  was  close 
to  u . .   Their  model  applied  to  the  mechanically  smooth  dense 
phase  fluidization  although  they  did  observe  bubble  formation 
and  slugging  in  the  dense  phase  region  of  fluidization. 

The  model's  motion,  Figure  16,  is  primarily  upward 
at  the  center  of  the  bed  and  downward  along  the  wall . 
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Temperature  measurements  showed  the  entire  bed  was  at  a  con- 
stant temperature  except  for  a  small  region  at  the  bottom 
and  a  narrow  layer  near  the  wall.   Heat  is  tranf erred  radially 
through  a  thin  air  film  to  the  particles  and  air  moving  down- 
ward outside  this  film.   The  particles  carry  the  heat  to  the 
bottom  of  the  bed  where  thermal  equilibrium  is  attained  almost 
instantaneously  with  the  incoming  cold  air.   The  bombardment 
of  small  particles  prevents  the  formation  of  a  laminar  bound- 
ary layer,  and  leaves  a  very  thin  laminar  sublayer  and  a 
thicker  turbulent  outer  layer. 

The  heat  transfer  rate  from  the  surface  to  the  fluid- 
ized  bed  was  found  to  be  independent  of  the  thermal  conduc- 
tivity of  the  solid  particles,  and  the  large  isothermal 
region  of  the  fluidized  bed  could  be  destroyed  if  the  down- 
ward flow  of  the  small  particles  was  short-circuited  into 
the  main  flow  before  they  reached  the  bottom,  i.e.  under 
slugging  conditions . 

Levenspiel  and  Walton  [41]  presented  a  similar  model 
where  the  resistance  to  the  heat  flow  is  due  to  the  laminar 
gas  layer  which  is  destroyed  by  solid  particles  passing 
through  it.   Thus  the  average  thickness  of  the  laminar  layer 
is  much  less  than  in  an  empty  tube.   They  assumed  that  the 
particles  are  of  uniform  diameter  (d  ) ,  stationary,  and 
arranged  in  equally  spaced  horizontal  layers.   The  boundary 
layer  formation  at  the  wall  is  similar  to  that  of  a  flat 
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plate,  see  Figure  17.   At  points  of  contact  with  the  station- 
ary particles,  the  boundary  layer  is  destroyed  and  starts 
once  again. 

Unfortunately,  the  two  models  do  not  take  into  account 
the  influence  of  the  solid  particles  on  the  heat  transferred, 
so  they  must  be  considered  of  limited  application. 

2 .   Unsteady  Heat  Conduction  by  Single  Particles  in 
Direct  Contact  with  the  Heat  Exchange  Surface 

With  the  advent  of  photographic  studies  of  the  bound- 
ary wall  during  f luidization,  it  was  shown  that  the  particles 
tended  to  associate  in  groups  which  moved  slowly  along  the 
wall  or  remained  in  contact  with  the  wall  a  certain  length  of 
time.   Those  observations  led  Botterill  et  al .  [6]  to  propose 
that  heat  was  transferred  from  the  vertical  surface  to  the 
adjacent  layer  of  particles  arriving  together  at  the  wall. 
If  the  residence  time  of  the  particle  near  the  wall  was  long, 
heat  would  penetrate  further  into  the  bulk  of  the  particles 
and  if  the  residence  time  was  short,  heat  exchange  calcula- 
tions would  be  limited  to  only  one  particle,  the  nearest  to 
the  surface.   They  confined  their  theoretical  considerations 
to  an  isolated  particle  lying  in  contact  with  the  wall. 

After  much  experimentation  this  model  was  extended 
to  cover  different  particle  sizes  and  materials  and  a  wide 
range  of  residence  times  for  the  particles  at  the  surface. 
The  results  showed  the  validity  of  the  model  for  short  resi- 
dence times  with  the  observed  effect  that  the  heat  transfer 
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coefficient  is  a  decreasing  function  of  the  residence  time. 

Again,  this  model  can  only  be  applied  to  mechanically  smooth 

dense  phase  fluidization  as  the  model's  validity  breaks  down 

in  the  lean  phase  or  when  bubbling/slugging  occurs. 

3 .   Unsteady  Transfer  of  Heat  to  "Packets"  of  Particles 
Which  are  Renewed  by  Violent  Disturbances  in  the 
Core  of  the  Fluidized  Bed 

In  an  attempt  to  present  a  mechanism  for  the  heat 
transfer  in  bubbling,  dense  phase  beds,  Mickley  and  Fairbanks 
[47]  proposed  a  model  of  a  small  group  of  particles  moving  as 
individual  packets.   They  based  their  model  on  the  fact  that 
the  voidage  of  the  dense  phase  bubbling  beds  was  found  to 
be  close  to  the  voidage  of  quiescent  beds. 

In  the  dense  phase  fluidized  bed,  each  particle  was 
expected  to  be  in  contact  with  several  neighbors  most  of  the 
time.   The  packets  are  not  considered  permanent;  they  have  a 
finite  life  expectancy;  and  their  voidage,  density,  heat  con- 
ductivity, and  heat  capacity  are  assumed  to  be  the  same  as 
those  of  the  quiescent  bed.   See  Figure  18. 

This  model  assumes  that  as  a  packet  at  bed  tempera- 
ture T,  comes  into  contact  with  a  flat  surface  at  a  tempera- 
ture T   that  unsteady  heat  transfer  will  start  upon  contact, 
w 

Further,  if  A   is  the  contact  area  between  the  packet  and  the 
m 

wall,  and  the  packet  is  homogenous,  then  the  heat  transferred 
after  some  time  t  will  be  given  by  the  equation 

cr.=h.TA   (  T   -  T,  )  ,  (3.9) 

^wi     wiL   m    w    b 
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where  the  local  instantaneous  heat  transfer  coefficient  is 
given  by  the  equation 


h  .T 
wiL 


k   p   c 
m  m  m 


TT  t 


0.5 


(3.10) 


Mickley  and  Fairbanks  [47]  also  presented  two  devia- 
tions from  their  original  model  which  covered  the  phenomenons 
of  slug  flow  and  side  mixing.   The  first,  on  slug  flow,  was 
for  fluidized  beds  with  very  low  gas  flow  rates  and  no  turbu- 
lent mixing.   In  this  model,  Figure  19,  the  particles  were 
to  move  downward  along  the  outside  wall  of  the  bed.   The  flow 
is  not  smooth  but  is  spasmodic  owing  to  the  intermittent 
passage  of  ascending  gas  bubbles  near  by  in  the  core  of  the 
bed.   However,  assuming  that  the  particles  moved  downward  at 
a  constant  velocity  (u  ) ,  the  age  of  all  the  packets  at  a 
distance  I    from  the  bed  top  was  t  =  i/u      and  the  resulting 


average  heat  transfer  coefficient  for  a  heater  of  length  L 
could  be  found  from  the  equation 


H 


h   =  2 
w 


k   p   c   u 
m  m  m   s 


7T  L 


H 


0.5 


(3.11) 


The  second  model  deviation,  on  side  mixing,  was  for 
beds  with  large  heat  transfer  surfaces  that  were  highly 
turbulent.   To  represent  this  case,  Figure  20,  they  assumed 
that  as  the  particle  packet  is  moving  downward  some  of  it 
is  exchanged  with  some  other  packet  at  temperature  T,  ,  brought 
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in  sideways  from  the  core  of  the  bed.   The  average  number  of 
replacement  packets  at  the  wall  per  unit  time  resulting  from 
side  mixing  was  then  defined  as  s,  resulting  in  the  local  heat 
transfer  coefficient  equation 

h  T  =  (  k   o   c   s  )  °-5  .  (3.12) 

wL      m  m  m 

It  should  be  noted  that  for  these  penetration  type 
models  the  residence  time  of  the  particle  is  of  fundamental 
importance  and  has  generated  a  considerable  amount  of  research 
both  in  its  determination  and  its  effect  on  the  various  models 
of  this  section. 

Patel  [53],  in  a  further  spin-off  of  Mickley  and 
Fairbank's  original  models  presented  two  surface  renewal 
models  that  can  be  characterized  as  follows. 

In  model  one,  Figure  21,  a  packet  initially  at  the 
bed  temperature  T,  arrives  at  the  wall  (temperature  T  )  at 

c  D  W 

time  t  =  0.   The  packet  is  assumed  to  have  the  same  proper- 
ties as  the  bed  at  minimum  f luidization.   It  receives  heat 
from  the  wall  through  a  contact  resistance  and  returns  to 
the  bed  after  some  length  of  time.   During  that  residence 
time,  heat  penetrates  a  distance  x;  beyond  this  distance 
the  temperature  is  constant  and  equal  to  the  bed  tempera- 
ture T,  . 
b 

In  model  two,  Figure  22,  the  fluidized  particles 

are  assumed  to  be  spherical  and  of  the  same  diameter  (d  ) . 

P 

A  particle  at  the  bed  temperature  T,  arrives  at  the  heating 
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Figure  21.   Pate  1  Model  I 


64 


surface  (temperature  T  )  and  receives  heat  from  the  fluid 

w 

adjacent  to  the  wall.   This  fluid  is  assumed  to  be  stagnant 

and  at  wall  temperature  T  .   While  the  particle  is  at  the 

w  *■ 

surface  it  also  loses  heat  by  conduction  to  a  packet  of  par- 
ticles of  thickness  x  situated  between  the  wall  particle 
and  the  bulk  of  the  bed.   The  packet  is  assumed  to  have  the 
same  properties  as  the  bed  at  minimum  fluidization  conditions. 

Patel  by  no  means  had  the  only  variation  on  Mickley 
and  Fairbank's  packet  theory.   Numerous  other  variations 

exist  with  different  h   values  being  correlated  by  each 

w  3  2 

researcher.   Appendices  A  through  C  give  some  of  the  result- 
ing correlations. 

4 .   Steady  Convection  Through  the  Emulsion  Layer 

van  Heerden  et  al .  [66]  observed  that  the  heat  capa- 
city of  the  particles  per  unit  volume  is  of  an  order  of  mag- 
nitude of  thousands  times  greater  than  that  of  the  gas  and 
that  the  mean  particle  velocity  is  much  lower  than  the  gas 
velocity.   Therefore,  the  largest  portion  of  the  heat  will  be 
transferred  by  the  moving  particles.   They  assumed  that  the 
gas  temperature  was  that  of  the  particles  and  as  a  result 
that  the  gas  only  provided  a  suspension  for  the  particles. 

They  observed  negligible  radial  motion  of  the  bed 
particles  and  therefore  concluded  that  the  radial  heat  trans- 
fer was  a  function  of  the  thermal  conductivity  of  the  sus- 
pension.  Experimentally,  they  verified  that  the  heat  transfer 
coefficient  is  large  for  short  heat  exchange  surfaces  and  vice 

versa. 
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Wicke  and  Fetting  [75]  furthered  van  Heerden ' s  work 

and  proposed  a  model,  Figure  23,  for  those  observations.   In 

their  model,  heat  (q  )  from  an  exchange  surface  was  first 

w 

transferred  by  conduction  through  a  gas  layer  whose  thickness 

was  5-.   This  heat  was  then  divided  into  two  components: 

namely  heat  (q  )  taken  by  the  particles  flowing  parallel  to 

the  surface  in  a  second  zone  of  emulsion  of  thickness  5  ,  and 

e 

heat  (q  )  transferred  into  the  core  portion  of  the  bed  by 
interchange  of  the  particles.   In  their  calculations,  how- 
ever, they  neglected  the  term  q  . 

The  last  model  to  be  presented  is  actually  a  combina- 
tion of  several  heat  transfer  mechanisms  proposed  in  earlier 
models . 

5 .   Kunii  and  Levenspiel  Model  [34] 

In  this  model,  Figure  24,  four  mechanisms  operate 
together . 

a.  A  film  of  gas  coats  the  surface.  Its  thickness 
(<$r)  is  large  or  small  depending  on  whether  a  bubble  of  air 
is  near  the  surface  or  the  emulsion  is  uniform  and  close  to 
the  surface. 

b.  Some  solid  particles  are  in  direct  contact  with 
the  heat  exchange  surface. 

c.  There  is  a  layer  of  emulsion  with  thickness  5 

1  e 

which  flows  along  the  wall. 

d.  Part  of  the  emulsion  layer  is  replaced  occasion- 
ally by  fresh  emulsion  coming  from  the  core  of  the  bed  or  by 

bubbles  rising  along  the  wall. 
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Figure  23.   Wicke  and  Fetting  Model 
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Depending  upon  the  fluidization  conditions  and  the 
position  of  the  heat  exchange  surface  in  the  bed,  emphasis 
can  be  placed  upon  any  combination  of  the  above  four 
mechanisms . 

C.   FLUIDIZED  BED  VARIABLES  THAT  AFFECT  HEAT  TRANSFER 

In  addition  to  the  large  number  of  models  that  have  been 
proposed  by  various  researchers  to  explain  their  experimental 
results  in  heat  transfer,  there  are  a  number  of  variables 
which  have  been  found  to  affect  the  heat  transfer  rate.   In 
general,  the  factors  which  will  control  the  heat  transfer 
characteristics  of  the  bed  will  be  those  variables  which 
determine  the  fluidizing  properties  of  the  bed.   In  this 
section  an  attempt  will  be  made  to  list  some  of  the  more 
important  factors  and  their  effect  on  heat  transfer.   An 
excellent  indepth  review  of  general  fluidized  bed  character- 
istics was  researched  by  Gutfinger  and  Abauf  [28] .   Addition- 
ally, Appendices  A  through  C  present  a  list  of  nondimensional 
correlations  reported  by  different  investigators  while  trying 
to  correlate  their  experimental  findings. 
1.   Properties  of  the  Fluid 

a.   Density  (pf)  and  Specific  Heat  (c  f)  of  the  Fluid 
The  dependence  of  the  heat  transfer  coefficient 
on  the  specific  heat  of  the  gas  has  not  been  clearly  deter- 
mined as  the  existing  data  shows  no  common  trend.   However, 
Gelperin  and  Einstein  [18]  have  predicted  that  there  would 
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be  an  increase  in  the  heat  transfer  rate  with  an  increase  in 
c  -  or  the  product  pfc  -  when  the  fluid  was  at  high  pressures 
and  velocities. 

b.  Gas  Viscosity  (u) 

The  majority  of  researchers  agree  that  with  an 
increase  in  the  gas  viscosity  the  wall  heat  transfer  coeffi- 
cient will  decrease. 

c.  Thermal  Conductivity  of  the  Gas  (kf) 

This  is  the  fluid  property  that  has  the  most  pro- 
nounced effect  on  the  heat  transfer  coefficient.   Depending 

on  which  researcher  is  consulted  the  value  of  h  will  increase 

w 

proportional  to  kf  to  the  power  of  0.5  -  0.67. 
2 .   Properties  of  the  Bed  Particles 

a.  Diameter  (D  )  and  Shape  of  the  Particle  (6) 

P  s 

A  significant  number  of  researchers  [37,47]  have 
concluded  that  the  heat  transfer  coefficient  decreases  when 
large  diameter  particles  are  used.   However,  the  results  of 
Sarkit  et  al .  [60]  as  reported  by  Gelperin  and  Einstein  [18] 
show  that  for  the  laminar  flow  region  the  heat  transfer 
coefficient  varies  inversely  with  particle  diameter  and  for 
the  turbulent  flow  region  it  has  a  direct  variation. 

The  values  for  the  heat  transfer  coefficient  tend 
to  increase  as  the  sphericity  (<j>  )  increases  (rounder  and 
smoother  particles)  . 

b.  Density  of  the  Particle  (p  ) 

Experimentation  has  shown  that  the  wall  heat 

transfer  coefficient  increases  with  the  increasing  density 
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of  the  solid  particles.   However,  this  increase  slows  in 
its  magnitude  as  the  fluid  flow  becomes  more  turbulent. 

c.  Specific  Heat  of  the  Particle  (c   ) 

r  ps 

The  heat  transfer  coefficient  increases  with  an 
increase  in  the  specific  heat  of  the  particle. 

d.  Thermal  Conductivity  of  the  Particle  (k  ) 
Based  on  measured  experimental  results  over  a 

wide  range  of  particle  thermal  conductivities  it  has  been 
determined  that  the  effect  of  the  thermal  conductivity  of 
the  particle  on  the  heat  transfer  coefficient  is  slight,  so 
much  so  that  it  is  usually  neglected  as  a  factor  of  concern. 
3 .   Fluidization  Conditions 
a.   Gas  Mass  Velocity  (G) 

The  increase  in  the  heat  transfer  coefficient  is 
proportional  to  G  where  n  will  vary  according  to  the  experi- 
ment being  conducted.   Its  range  has  been  shown  to  be  quite 
large.   About  the  only  facts  that  have  been  agreed  upon  by 
researchers  for  this  variable  is  that  of  the  general  shape 
of  the  gas  mass  velocity  versus  heat  transfer  coefficient 
curve  [72],  Figure  25,  and  that  there  is  some  maximum  value 

for  h,  which  exists  at  some  gas  mass  velocity  G   ..   Tables 

w  3  2      opt 

listing  the  nondimensional  correlations  for  both  the  rising 

and  falling  sides  of  the  curve,  for  the  maximum  heat  transfer 

coefficients,  and  for  the  special  gas  flow  rates  when  the 

h ,  „  is  attained  are  presented  by  Einstein  and  Gelperin 

wmax  r  j.  r- 

[14]. 
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Figure  25.   Heat  transfer  coefficient  vs.  mass  velocity. 
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b.   Bed  Porosity  or  Voidage  (e) 

Based. on  film  theory  the  heat  transfer  coeffi- 
cient will  vary  inversely  with  some  function  of  the  bed 
porosity.   That  is,  as  the  particle  aggregate  at  the  wall 
decreases  so  will  the  heat  transfer  rate.   However,  there 
has  been  some  experimental  data  which  has  demonstrated  the 
opposite  effect. 

4 .   Geometric  Variables 

a.  Bed  Diameter  (D) 

It  has  not  been  demonstrated  that  there  exists 
a  dependency  of  the  heat  transfer  coefficient  on  the  bed 
diameter . 

b.  Diameter  of  an  Immersed  Heater  (D„) 

n 

Researchers  are  in  general  agreement  that  the 
wall  heat  transfer  coefficient  increases  with  a  decrease  in 
the  immersed  heater  (tube)  diameter,  and  that  the  heat  trans- 
fer rate  diminishes  as  the  diameter  of  the  heater  is  increased 
These  effects  are  valid  up  to  a  particle  diameter  of  10  mm 

(0.39  in)  after  which  h,  is  independent  of  the  heater 

w        r 

diameter . 

c.  Fluidized  Bed  Height  (L) 

There  has  been  no  dependence  shown  between  the 
heat  transfer  coefficient  and  the  fluidized  bed  height. 

d.  Length  of  Heat  Exchange  Surface  (L„) 

This  is  a  variable  that  has  had  several  differ- 
ent results  on  its  effect  on  the  heat  transfer  coefficient 
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all  of  which  have  been  substantiated  by  additional  research. 
Basically,  one  research  group  has  found  that  the  heat  trans- 
fer coefficient  decreases  with  an  increase  in  the  length  of 
the  heat  exchange  surface.   Another  research  group  concluded 
that  if  LH/D  was  greater  than  7  for  an  external  heat  transfer 
surface,  then  the  heat  transfer  coefficient  was  independent 
of  heater  length  and  that  internal  heat  exchange  surfaces  had 
no  observed  effect  on  the  heat  transfer  coefficient.   And  yet 
a  third  research  group  found  that  if  the  heat  exchange  surface 
was  located  above  the  active  section  of  the  bed  (mixing 
region) ,  then  the  heat  transfer  coefficient  would  be  inde- 
pendent of  the  heater's  length. 

e.   Heat  Exchange  Surface  Roughness 

It  has  been  determined  [31,68]  that  the  heat  trans- 
fer coefficient  is  sensitive  not  only  to  the  type  and  orient- 
ation of  the  roughness  but  also  to  the  size  of  the  particles 
relative  to  the  size  of  the  roughness.   When  the  particles 
are  of  a  large  size  relative  to  the  roughness,  the  gas  con- 
duction paths  will  be  longer  and  the  heat  transfer  will  be 
lower  as  compared  to  a  smooth  tube  or  surface.   See  Figure 
26a  and  26b. 

When  particle  and  roughness  sizes  are  similar, 
the  multiple  contact  should  reduce  the  mean  gas  conduction 
paths  and  gives  rise  to  a  considerable  increase  in  the  heat 
transfer,  Figure  26c.   However,  possibly  due  to  particle 
capture   and  the  consequent  increase  in  particle  residence 

74 


smooth    sur-face 


a 


fine  roughness 

Q 


(a) 


(b) 


matching  roughness 


coarse  roughness 


Figure  26.   Contact  geometry  of  surface-particle 
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times,  the  increase  in  the  maximum  heat  transfer  coefficient 
for  vertical  grooves,  is  of  the  order  of  50%  for  the  smaller 
particles  and  18%  for  the  larger  particles.   The  particle 
capture  and  inability  of  the  bubbles  to  replenish  particles 
in  the  case  of  horizontal  grooves  undoubtedly  contributes  to 
this  difference. 

Where  the  size  of  the  particles  is  much  smaller 
relative  to  the  roughness  size,  particle  capture  is  not  a 
serious  problem,  but  there  is  some  improvement  in  particle- 
surface  contact  as  well  as  a  larger  effective  area  for  heat 
transfer,  Figure  26d. 

5 .   Miscellaneous  Variables 

a.  Entrance  Effects 

Entrance  effects  in  a  flow  situation  cannot  be 
ignored  or  avoided.   Control  or  limiting  entrance  effects  can 
only  be  accomplished  by  careful  design  of  the  gas  distributor, 
which  in  turn,  will  influence  the  configuration  of  the  active 
(mixing)  region  of  the  fluidized  bed. 

b.  Temperature  and  Pressure  of  the  Fluidized  Bed 
The  heat  transfer  coefficient  is  independent  of 

the  pressure  of  the  fluidized  bed  or  fluidizing  fluid.   The 
temperature  of  the  fluid,  however,  is  another  matter.   In 
general  it  can  be  assumed  that  the  heat  transfer  coefficient 
is  independent  of  the  temperature  of  the  bed  or  fluid. 
Studies  by  Zabrodsky  [79]  and  Yoshida  et  al.  [77],  however, 
indicate  that  the  radiation  component  of  the  overall  heat 
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transfer  coefficient  must  be  taken  into  account  when  the  bed 
temperature  exceeds  900-1300  K.   In  turn,  the  magnitude  of 
the  radiation  component  effect  will  be  dependent  on  the  bed 
container  geometry. 

D.   IMMERSED  BODIES 

Experimental  work  conducted  in  fluidized  beds  can  be  sub- 
divided into  two  distinct  categories.   The  first  deals  with 
all  experimental  work  that  uses  an  external  heat  transfer 
surface,  i.e.  the  wall  of  the  fluidized  bed  container  itself 
is  used  as  the  heat  exchange  surface.   The  second  category 
uses  an  internal  heat  transfer  surface.   In  this  case,  bodies 
of  different  shapes  (sphere,  cylinder,  plate,  coil,  etc.) 
which  are  immersed  in  the  bulk  of  the  fluidized  bed  become 
the  heat  transfer  surface. 

While  the  previously  discussed  variables  and  character- 
istic bed  behavior  are  equally  germane  to  either  category, 
the  use  of  immersed  heating  surfaces  will  have  additional 
factors  that  will  need  to  be  taken  into  account  during  a 
design  process.   The  location  of  the  heat  source  with  respect 
to  the  gas  distributor  plate  and  its  plane  of  orientation 
will  have  a  significant  effect  on  the  heat  transfer  coeffi- 
cient.  The  results  of  the  variations  are  sometimes  contra- 
dictory and  depend  strongly  on  the  local  hydrodynamic 
conditions  in  each  case. 

Limited  experimentation  has  been  conducted  in  the  use  of 
spheres,  flat  plates  and  small  cylinders  with  the  data 
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obtained  being  sketchy  at  best.   The  bulk  of  the  immersed 
body  experimentation  has  been  with  tubes  for  which  additional 
variations  in  heat  transfer  tube  design  (size,  shape,  spacing, 
gap,  pitch  and  material)  will  complicate  the  general  heat 
transfer  mechanism.   An  excellent  indepth  coverage  of  this 
topic  was  conducted  by  Saxena  et  al .  [61]  and  is  summarized 
as  follows. 

1 .   Horizontal  Tubes 

The  total  heat  transfer  coefficient  (h  )  for  single 

w         3 

cylindrical  tubes  increases  with  factors  that  increase  the 

gas  thermal  conductivity  near  the  tube  wall  (i.e.  temperature 

at  the  tube  surface)  and  reduces  the  particle  residence  time 

(i.e.  decrease  in  tube  diameter  and  increase  in  fluidization 

velocity) .   However,  at  higher  gas  velocities  the  decrease 

in  bed  density  reduces  the  heat  absorption  capacity  of  the 

fluidized  particles,  resulting  in  a  decrease  in  the  heat 

transfer  coefficient  (h  ) .   The  heat  transfer  coefficient 

w 

also  increases  with  decreasing  particle  diameter  (d  ) . 

The  heat  transfer  coefficient  is  highest  for  tubes 
with  a  lenticular  cross  section  for  dense  moving  packed  beds. 
Gas  pockets  are  not  created  in  the  flow  wake  (downstream 
side)  of  a  lenticular  tube,  and  the  particles  are  not  stag- 
nant on  the  upstream  side  of  the  tube. 

The  presence  of  surrounding  tubes  in  an  array  affects 
the  total  heat  transfer  coefficient  for  an  individual  tube 
by  resisting  particle  movement  and  by  altering  the  particle 
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temperature  from  that  of  the  bed  bulk,  which  reduces  the  tem- 
perature driving  force.   This  is  at  least  one  of  the  reasons 
why  the  heat  transfer  coefficients  for  the  lower  tubes  are 
greater  than  for  the  upper  tubes,  and  a  close  spacing  of  the 
tubes  results  in  a  decrease  of  the  heat  transfer  coefficient. 

Fins  will  increase  the  rate  of  heat  transfer  to  the 
tubes.   However,  their  effectiveness  increases,  reaches  a 
maximum,  and  then  decreases  with  decreasing  spacing  distance 
and  increasing  fin  length  [54]. 

An  increase  in  tube-to-tube  spacing  in  a  bundle  of 
short-finned  tubes,  which  reduces  the  resistance  to  particle 
movement  between  the  tubes,  increases  the  total  heat  transfer 
coefficient.   However,  an  increase  in  tube-to-tube  spacing 
for  a  bundle  of  long-finned  tubes  does  not  noticeably  in- 
crease the  heat  transfer  rate,  implying  that  the  resistance 
to  particle  movement  between  fins  predominates. 
2  .   Slanted  Tubes 

There  does  not  appear  to  be  any  particular  advantage 
in  slanting  the  heat  transfer  tubes.   Minimum  heat  transfer 
coefficients  have  been  observed  at  45°  for  bare  tubes  and  at 
60°  for  finned  tubes  with  respect  to  the  horizontal  bed  axis. 
The  heat  transfer  coefficient  values  obtained  were  in  compari- 
son to  the  values  obtained  for  tubes  of  the  same  configura- 
tion but  in  a  horizontal  and  vertical  position  [23] . 
3 .   Vertical  Tubes 

The  total  heat  transfer  coefficient  for  a  single  ver- 
tical tube  generally  increases  with  a  decrease  in  tube  length 
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and  diameter,  a  decrease  in  the  fluidized  particle  size  (sizes 
less  than  about  1  mm  (0.039  in)  at  atmospheric  pressure) ,  and 
an  increase  in  fluidization  gas  velocity  up  to  u.  . 

The  heat  transfer  coefficient  for  a  tube  within  a 
bundle  increases  with  an  increase  in  the  tube  spacing  and 
is  greater  for  tubes  located  nearer  to  the  center  of  the 
bundle,  at  least  for  small  bundles. 

The  heat  transfer  coefficient  for  a  finned  tube 
generally  increases  with  the  fin  height  and  a  decrease  in 
the  frequency  of  fins  along  the  tube  length  (in  the  range  of 
21-10  fins/in) .   Heat  transfer  also  improves  with  roughened 
and  grooved  tubes . 
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IV.   EXPERIMENTAL  APPARATUS  AND  PROCEDURES 

This  section  will  describe  the  details  and  procedures  as 
related  to  the  experimental  apparatus.   The  methodology  and 
experimental  procedures  will  also  be  discussed. 

A.   DESCRIPTION  OF  THE  FLUIDIZATION  APPARATUS 

The  fluidization  container  for  the  experimental  investi- 
gation was  constructed  of  one  half  inch  thick  clear  plexiglas 
(Grade  GM) .   The  inside  dimensions  of  the  rectangular  enclo- 
sure were  variable  as  one  wall  was  capable  of  movement  within 
the  boundaries  of  the  container.   This  movement  allowed  the 
variation  of  the  inside  dimensions  from  six  inches  wide  by 
six  inches  deep  by  eighteen  inches  high  to  twenty  four  inches 
wide  by  six  inches  deep  by  eighteen  inches  high  above  the 
distributor  plate.   All  joints  were  glued  together  and  wood 
screws  added  in  mounting  the  distributor  plate  for  additional 
weight  carrying  support.   Mounted  within  each  of  the  two  side 
walls  is  a  five  inche  wide  by  ten  inch   high  strip  heater 
that  has  the  heating  surface  on  the  same  plane  as  the  inside 
plexiglas  wall.   Backing  the  strip  heaters  are  two  layers  of 
insulation  (silicone  rubber  sponge  and  plexiglas)  which 
provide  both  rigidity  and  a  mounting  platform  for  the  strip 
heater.   Static  pressure  taps  and  bed  thermocouple  taps  were 
installed  at  evenly  spaced  intervals  as  shown  in  Figure  27. 
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Figure  27.   Fluidization  apparatus 
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Pressure  drops  were  measured  using  the  static  pressure  taps 
with  the  output  going  to  a  Meriam  Type  W  Model  A324  sixty 
inch  water  manometer  or  a  Model  33KA35  thirty  inch  water 
manometer  via  a  pressure  manifold.   Air  and  bed  temperatures 
within  the  experimental  apparatus  are  measured  with  copper- 
constantan  thermocouples  as  described  in  the  EXPERIMENTAL 
HEAT  TRANSFER  PROCEDURE  AND  APPARATUS  section  of  this  thesis 
and  recorded  on  a  Newport  Model  267B-TF2,  Multipoint  (12 
Channel)  Digital  Temperature  Recorder.   The  complete  equip- 
ment listing  for  this  experimental  study  is  given  in  Appendix 
E. 

The  distributor  plate  utilized  in  the  fluidization  con- 
tainer was  constructed  of  one  half  inch  thick  plexiglas 
(Grade  GM)  with  a  grid  of  0.125  in  diameter  holes  bored  at 
0.25  in  center  to  center  spacing.   To  cause  the  required 
pressure  drop  which  aids  the  air  distribution  and  to  prevent 
the  particles  from  falling  into  the  air  box,  one  layer  of  140 
mesh  stainless  steel  wire  cloth  was  mounted  on  top  of  the 
distributor  plate. 

The  air  was  supplied  by  a  Clements  Cadillac  Blower/ 
Suction  Cleaner  Model  G12  with  a  nominal  output  of  about  20 
scfm  of  air  at  a  nominal  line  pressure  of  1.5  psig.   The  flow 
output  data  is  based  on  experimental  testing  of  the  blower 
since  blower  characteristic  curves  are  not  available.   The 
air  inlets  to  a  two  and  one  half  inch  inside  diameter  plexi- 
glas tube  that  has  one  half  inch  circular  air  ports  spaced 
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at  one  inch  center  to  center  intervals  along  its  length.   The 
air  ports  are  located  at  the  underside  of  the  tube  allowing 
the  air  to  exhaust  downward  into  the  air  box.   The  plexiglas 
tube  runs  the  length  of  the  air  box  which  has  inside  dimen- 
sions of  twenty  four  inches  wide  by  six  inches  deep  by  six 
inches  high.   The  air  is  unfiltered  and  prevented  from  exit- 
ing outside  of  the  variable  bed  inlet  area  by  an  air  baffle 
mounted  to  the  movable  wall  above  the  distributor  plate.   The 
schematic  line  diagram  of  the  air  supply  system  is  shown  in 
Figure  28. 

The  air  flow  rate  was  measured  with  a  Fischer  and  Porter 
Company  Model  10A3565A  Flowmeter  (rated  at  11.1  scfm)  which 
was  corrected  for  local  barometric  pressure.   Photographs 
showing  two  views  of  the  overall  experimental  apparatus  can 
be  found  in  Figure  29  and  30. 

B.   EXPERIMENTAL  HEAT  TRANSFER  PROCEDURE  AND  APPARATUS 

Four  different  rectangular  bed  configurations  using  a 
standard  static  bed  height  were  investigated.   For  each  bed 
configuration  both  an  unheated  and  heated  run  were  conducted 
with  various  pressure,  temperature  and  visual  data  taken  on 
each  run. 

Once  the  bed  size  was  adjusted,  the  bed  of  particles  was 
initially  fluidized  at  a  point  above  minimum  f luidization . 
When  steady  state  conditions  had  been  reached,  the  bed  and 
air  (inlet  and  outlet)  temperatures  were  recorded  and  a 
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Figure    28.       Diagram    o-f    f 1 u i d i zat i on    system 
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Figure  29.   Photograph  of  experimental  apparatus 


86 


Figure  30.   Photograph  o-f  experimental  apparatus. 


87 


comment  recorded  describing  the  visual  condition  of  the  exter- 
nal boundaries  of  the  bed.   This  was  the  starting  point  for 
the  run. 

As  a  separate  study,  bed  pressure  drop  was  recorded  as  a 
function  of  several  flow  rates.   This  information  was  used  to 
determine  the  point  of  minimum  fluidization  for  the  particular 
bed  configuration  and  to  serve  as  an  indicator  for  the  pre- 
sence of  slugging  or  channeling  conditions  within  the  bed. 

After  the  base  data  run  was  recorded,  power  was  applied 
to  the  two  Watlow  Silicone  Rubber  strip  heaters  by  way  of  a 
Lamba  Regulated  Power  Supply,  Model  LK34  5A.   The  power  supplied 
to  the  heaters  was  computed  by  measuring  the  voltage  and 
amperage.   Energy  losses  by  conduction  were  accounted  for  in 
the  heater  configuration  and  an  energy  balance  computed  for 
each  run. 

Since  the  purpose  of  this  experiment  did  not  include  a 
need  to  determine  a  heat  transfer  coefficient,  provisions 
were  not  made  to  determine  heater  wall  temperatures  or  the 
temperature  gradients  close  to  the  heater  wall.   On  all  runs 
the  expanded  bed  height  covered  the  strip  heaters  ensuring 
that  all  heat  transfer  was  either  to  the  bed  or  losses 
through  the  insulation.   Radiation  losses  were  neglected 
since  bed/wall  temperatures  were  considerably  less  than  1000 
K. 

At  various  heater  power  levels  a  series  of  temperature 
readings  were  recorded  until  flow  characteristics  were 
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determined  for  each  bed  configuration.   The  above  experimental 
procedures  were  continued  until  all  bed  geometry  configura- 
tions had  been  studied. 

The  details  of  the  heater  configuration  and  mounting 
method  are  shown  in  Figure  31 .   Insulating  end  blocks  were 
used  for  both  mounting  and  end  heat  loss  reduction  for  the 
two  heaters.   Heater  insulation  consisted  of  0.5  in  of  sili- 
cone rubber  sponge  with  a  mounting  block  of  0.5  in  of  clear 
(Grade  GM)  plexiglas .   Two  insulated  copper-constantan  thermo- 
couples were  applied  to  the  outside  of  the  plexiglas  mounting 
blocks  to  measure  their  outer  surface  temperature.   These 
temperatures  were  used  to  estimate  the  heat  loss  through  the 
back  of  the  heater  by  conduction  for  use  in  the  heat  balance 
of  the  model. 

The  strip  heater  resistance  element  is  constructed  of  a 
fine  nicket  alloy  wire  wrapped  around  a  glass  cord.   The 
element  wires  are  arranged  in  a  parallel  rectangular  configur- 
ation thus  confining  any  failure  to  a  single  element  rather 
than  the  entire  unit.   Each  heater  is  rated  at  250  watts, 
115  volts  with  a  temperature  range  of  0-200°F.   A  consider- 
ably higher  temperature  range  is  possible  by  cementing  the 
heater  to  various  metal  plates,  but  is  limited  in  this  applica- 
tion due  to  the  melting  temperature  of  the  plexiglas. 

To  measure  the  bed  temperatures  a  grid  pattern  of  thirty 
retractable  Omega  Engineering  ANSI  T  (Copper-constantan) 
twenty-four  inch  stainless  steel  sheathed  exposed  junction 
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thermocouple   probes  was  used.   By  varying  the  depth  of  the 
individual  thermocouples  the  boundary  positions  of  the  par- 
ticle and  fluid  flow  cells  could  be  determined  and  changes 
to  the  original  pattern  noted  during  the  various  runs. 

C.   EXPERIMENTAL  APPARATUS  CALIBRATION 

Manometers  used  for  static  pressure  drop  and  differen- 
tial pressure  measurements  were  continuously  rechecked  for 
zero  readings  and  readjusted  as  required. 

The  temperature  recorder  and  thermocouples  were  calibra- 
ted with  a  Rosemount  Commutating  Bridge  Model  920A  and 
Rosemount  Calibration  Bath  Model  913A  using  liquid  nitrogen 
and  water.   A  least  squares  linear  curve  fit  was  applied  to 
the  resulting  data  for  each  thermocouple  to  obtain  suitable 
constants  for  temperature  correction. 

Technical  data  on  the  bed  particles  and  experimental 
verification  of  the  manufacturer's  specifications  can  be 
found  in  Appendix  F. 
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V.   PRESENTATION  AND  DISCUSSION  OF  RESULTS 

Based  on  bed  observations  during  the  fluidizing  process 
a  very  definite  electrostatic  charge  built  up  within  the  bed 
container,  see  Figure  22.   This  phenomenon  was  characterized 
by  the  bed  particles  adhering  to  the  walls  of  the  bed  after 
the  air  had  been  turned  off  and  to  the  thermocouple  probes 
during  f luidization.   Based  on  evidence  provided  by  Miller 
and  Logwinuk  [49]  this  phenomenon  will  to  some  extent  con- 
tribute to  erratic  results  in  the  heat  transfer  and  data 
measurements.   Some  of  the  factors  which  are  thought  to 
affect  the  magnitude  of  the  electrostatic  charge  effects 
are  (1)  the  dielectric  constant  of  the  bed  material  used, 
(2)  the  contact  potential  difference,  (3)  the  ratio  of  fluid- 
ized  bed  particles  to  gas,  (4)  the  properties  of  the  gas,  (5) 
the  temperature  of  the  solid-fluid  system,  and  (6)  the  nature 
of  the  surface  of  the  bed  particles.   For  the  purpose  of  this 
thesis,  the  electrostatic  charge  existance  will  be  acknow- 
ledged but  no  data  correction  will  be  considered  warranted. 
Further  study  on  this  area  is,  however,  considered  worthwhile 

A.   FLOW  PATTERNS 

The  particle  movement  along  the  wall  at  any  level  is  in- 
fluenced by  the  movement  within  the  fluidized  bed  at  that 
level,  as  it  is  the  movement  of  the  inner  particles  which 
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drag  the  outer  particles  along  the  wall.   Therefore  the  wall 
velocity  profile  should  represent  approximately  the  average 
particle  velocity  at  the  various  levels  within  the  bed.   Had 
the  motion  been  smooth  with  all  particles  starting  at  the 
top  of  the  bed  and  flowing  straight  down  to  the  bottom  of  the 
bed  then  the  ideal,  Figure  2,  flow  pattern  would  have  occurred 
and  the  velocity  of  the  particles  could  be  assumed  to  be  con- 
stant at  the  wall  with  a  positive  linear  gradient  into  the 
axis  of  symmetry  for  the  bed.   In  actuality  this  was  neither 
observed  nor  anticipated  within  the  rectangular  bed  under 
study  due  to  the  geometry  involved. 

In  the  initial  geometric  configuration,  that  of  a  square 
column  (6x6x12.5  in),  the  particle  motion  appeared  to  conform 
to  the  particle  circulation  pattern  of  Figure  33.   The  fluid- 
ization  somewhat  resembled  waves  in  that  there  were  intermit- 
tent fronts  moving  smoothly  up  through  the  bed  which  terminated 
in  a  percolating  (rupture)  action  at  the  top  of  the  bed.   The 
domed  bulge,  Figure  34,  of  particles  was  observed  roughly  in 
the  center  of  the  column  and  after  rupturing,  the  bed  particles 
would  be  thrown  to  the  sides  where  they  would  start  down  the 
bed  walls. 

It  was  not  possible  to  observe  any  distinct  points  of 
emergence  into  the  bed  of  the  particles  in  the  downward  tran- 
sit but  a  definite  entry  point  was  visible  at  the  bottom  of 
the  bed.   The  characteristic  shape  of  the  non-movement  area  was 
elliptical  along  the  bed  bottom  and  curving  up  at  the  corners , 

see  Figure  35. 
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Figure  32.   Electrostatic  charge 
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Figure  33.   Proposed  particle  cell  circulation. 
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Figure  34.   Square  -Muidized  column  -  top. 
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Figure  35.   Fluidized  column  -  bottom 
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Initially,  the  shape  of  the  non-movement  area  was  not 
elliptical.   If  the  side  of  the  apparatus  as  shown  in  Figure 
27  is  considered  to  be  the  front  face  then  bubbles  were  ini- 
tializing in  the  front  left  corner  of  the  bed  along  the 
movable  wall.   This  shift  in  symmetry  of  the  fluidizing  cell 
generation  was  resulting  in  a  shallower  non-movement  area  at 
the  left  front  corner  and  a  larger  build  up  at  the  right 
front  corner.   The  cause  of  the  problem  turned  out  to  be 
that  the  distributor  plate  was  not  at  a  right  angle  to  the 
air  flow.   Downward  motion  was  observed  at  all  four  walls  and 
the  non-movement  zone  was  located  within  one  inch  of  the 
bottom  of  the  bed. 

It  was  readily  observable  that  the  velocity  of  the  par- 
ticles emerging  at  the  center  of  the  bed  was  higher  than  that 
transiting  down  the  face  of  the  wall.   This  differential 
while  not  empirically  measured  during  this  study  was  attri- 
buted to  frictional  drag  at  the  wall  both  mechanical  and  as 
a  result  of  the  electrostatic  charge  present  at  the  bed  wall. 
There  is  some  justification  provided  by  Toomey  and  Johnstone 
[62].   Marsheck  [43]  also  gives  an  excellent  treatise  on  the 
velocity  profile  across  a  gas-f luidized  bed. 

The  height  of  the  domed  bulge  would  pulse  during  the 
fluidization  more  so  as  the  air  velocity  was  increased.   If 
the  air  blower  was  fully  energized  with  the  initial  bed  con- 
figuration then  the  entire  top  of  the  bed  would  rise  up  as 
a  mass  and  then  collapse  on  itself.   This  was  an  example  of 
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slug  flow  where  the  rising  bubbles  had  coalesced  into  one  large 
bubble  stretching  from  one  wall  to  another.   Slug  flow  was  not 
achieved  or  expected  in  any  of  the  other  configurations  as  the 
bed  was  not  deep  enough  nor  was  the  air  velocity  high  enough 
to  provide  the  right  conditions  for  the  bed  cross-sectional 
areas  present. 

For  the  air  flow  that  just  permitted  particle  circulation 
the  bed  height  fluctuations  were  negligible  although  the  par- 
ticles at  the  wall  still  moved  downward  in  an  irregular 
slightly  surging  movement.   It  is  felt  that  at  some  point 
in  the  bed  for  this  minimum  fluidization  velocity  that  the 
flow  pattern  internally  is  such  that  the  flow  pattern  dis- 
tributes the  particle  velocities  more  evenly  throughout  the 
column.   In  other  words,  the  circulating  motion  caused  by 
the  higher  velocity  in  the  center  and  slower  velocity  along 
the  wall  gives  way  to  a  more  random  or  eddying  motion  that 
evenly  distributes  the  small  air  bubbles  within  the  two 
velocity  boundaries  thus  preventing  rapid  bubble  coalescence. 
These  eddy  cells  once  formed  in  the  mixing  region  of  the  bed 
will  be  moved  up  through  the  bed  by  the  surges  resulting 
from  the  cell  rupture  at  the  bed  top.   Between  each  surge 
line  it  is  theorized  that  only  one  eddy  cell  will  exist  and 
that  it  will  grow  as  it  progresses  through  the  bed  until  it 
coalesces  with  the  gas  in  the  central  cell  core  at  the  bed 
top  prior  to  cell  rupture. 


97 


It  was  also  noted  that  as  the  thermocouple  probes  were 
lowered  into  the  bed  that  a  flattening  of  the  bed  top  would 
occur.   This  effect  became  more  pronounced  as  the  probe  pene- 
tration increased  and  is  felt  to  be  as  a  result  of  the  probe 
providing  a  grounding  channel  for  bubble/particle  flow.   In 
effect  the  probes  acted  as  a  calming  agent  in  the  agitation 
of  the  bed  without  changing  the  overall  bed  flow  cell 
configuration . 

When  the  movable  wall  was  positioned  for  a  6.25  inch  width 

(A  =37.5  in2)  and  also  for  an  8  inch  width  (A   =  48  in2)  the 
c  c 

bed  expansion  observed  at  the  minimum  fluidizing  air  velocity 
was  0.125-0.375  in.   When  the  bed  width  was  then  increased 
to  10  inches  (A  =  60  in2)  and  subsequently  12  inches  (A   = 
72  in2)  the  bed  expansion  observed  did  not  exceed  0.25  in. 
This  measurement  does  not  take  into  account  the  height  in- 
crease of  the  cell  dome  prior  to  rupture.   There  was  no 
general  horizontal  movement  of  the  particles  near  the  wall 
except  at  the  top  of  the  fluidized  bed  and  in  a  rather 
turbulent  mixing  region  near  the  bottom  of  the  fluidizing 
bed. 

As  the  movable  wall  was  repositioned  thus  increasing  the 
bed  air  inlet  area  the  bed  became  an  aggregative  fluidized 
bed.   This  resulted  in  multiple  cell  formations  within  the 
bed  vice  the  single  cell  observed  in  the  initial  bed  config- 
uration.    Unfortunately,  there  was  no  observable  pattern 
to  the  cell  formations  as  their  overall  positioning  was 
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Figure  36.   Bubbling  bed  photograph  (flc  =  60  in*). 


Figure  37.   Bubbling  bed  photograph  (Rc  =  72  in2) 
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random  in  nature,  see  Figures  36  and  37.   Only  one  cell  was 
observed  to  break  the  surface  of  the  bed  at  any  given  instant. 
Occasionally,  small  bubbles  were  noticed  breaking  the  surface 
of  the  bed  between  the  larger  bubble  cells.   This  was  thought 
to  be  caused  by  channeling  between  the  larger  bubble  cell 
boundaries . 

There  was  no  noticeable  change  in  the  outward  appearance 
of  the  particle  cell  flow  patterns  in  any  of  the  bed  configura- 
tions when  the  external  heaters  were  energized  during 
fluidization . 

B.   PRESSURE  DISTRIBUTION 

Static  pressure  readings  were  taken  at  various  intervals 
along  the  height  of  the  fluidized  bed  on  all  experimental  runs. 
In  order  to  determine  the  gas  mass  flow  rates  for  minimum 
fluidization,  measurements  of  the  pressure  drop  across  the 
bed  were  also  made  at  various  flow  rates.   All  data  obtained 
during  the  various  runs  can  be  found  in  Appendix  H .   The  mass 
flow  rate  at  the  maximum  pressure  drop  was  then  taken  as  the 
point  of  minimum  fluidization  in  each  case.   The  maximum 
pressure  drop  was  used  because  it  is  the  most  sharply  defined 
point  on  the  plots,  see  Figure  38. 

The  results  of  plotting  the  pressure  drop  versus  gas  mass 
velocity  profiles  produced  very  similar  patterns  which  are  in 
general  agreement  with  the  expected  curve  for  the  distributor 
design  used.   As  can  be  seen  from  Figure  38,  when  the  air  inlet 
bed  cross-sectional  area  was  increased  the  curve  would  shift 
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Figure  38.   Bed  pressure  drop  vs.  gas  mass  velocity. 
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to  the  left.   This  displacement  may  have  been  as  a  result  of 
the  lessening  in  the  contribution  of  the  mechanical  wall 
friction  drag  to  the  total  overall  drag  experienced  by  the 
bed  with  the  increased  air  flow  for  f luidization .   Neverthe- 
less, it  would  appear  as  the  bed  air  inlet  cross-sectional 
area  was  increased  the  gas  mass  velocity  required  for  fluid- 
ization  would  also  increase.   Additionally,  in  all  cases  ob- 
served, the  log  Ap  value  for  minimum  f luidization  was 
approximately  1.25. 

The  results  of  the  pressure  measurements  seemed  to  indicate 
that  mechanically  smooth  f luidization  was  taking  place. 
Channeling  and/or  slugging  was  not  occurring  at  the  minimum 
fluidization  operating  conditions  although  it  was  possible 
to  demonstrate  both  channeling  and  slugging  for  the  initial 
bed  configuration.   It  was  noted  that  during  fluidization 
small  fluctuations  were  visible  in  the  liquid  levels  in  the 
manometer  that  was  measuring  the  static  pressure  in  the  bed. 
This  operating  characteristic   is  in  agreement  with  similar 
conditions  observed  by  Parent  et  al .  [52]  in  their  study  of 
gas  fluidization. 

The  chain  of  events  in  measuring  the  pressure  drop  is  as 
follows.   Once  the  bed  was  initially  settled  the  air  flow 
was  initialized  to  the  bed.   As  the  air  flowed  through  the 
bed,  a  corresponding  pressure  drop  occurred.   This  pressure 
drop  increased  linearly  with  the  gas  flow  until  at  some  point 
the  bed  began  to  expand.   At  this  point  the  pressure  drop 
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became  nonlinear  yet  still  increased  to  some  maximum  value. 
This  maximum  value  was  picked  as  the  minimum  fluidization 
point  as  previously  stated.  The  pressure  drop  would  then 
decrease  until  the  bed  started  to  circulate  or  fluidize  at 
which  time  the  pressure  drop  would  stabilize  at  some  point 
lower  than  the  maximum  value  recorded  earlier.   The  pressure 
drop  remained  steady  as  the  air  velocity  was  increased  further 
but  soon  started  to  rise  again.   This  new  rise  in  the  pressure 
drop  across  the  bed  was  an  indication  that  slugging  was  start- 
ing to  set  in  and  was  especially  prevalent  in  the  initial  bed 
geometry.   It  was  also  noted  that  this  secondary  pressure 
drop  increase  was  less  prevalent  as  the  bed  cross-sectional 
area  was  increased. 

C.   TEMPERATURE  DISTRIBUTION 

Air  entered  the  air  box  at  temperatures  between  70°  and 
120 °F  and  in  every  case  left  the  bed  at  a  temperature  within 
1°F  of  that  of  the  top  of  the  bed,  as  measured  at  a  point  in 
the  axial  center  of  the  bed.   Only  two  data  runs  were  made 
with  the  bed  instrumented  to  record  the  temperature  distribu- 
tion.  The  first  run  (A   =37.5  in2)  had  an  average  air  inlet 
temperature  of  10 8° F  ±2.5°F  and  attained  a  constant  air  outlet 
temperature  of  125° F  ±0.4°F,  see  Figure  39.   The  second  tem- 
perature instrumented  data  run  (A  =  72  in2)  had  an  average 
air  inlet  temperature  of  116 °F  ±2°F  and  attained  a  constant 
air  outlet  temperature  of  129. 5°F  ±0.7°F,  see  Figures  40A  and 
40B. 
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In  the  graphics  of  Figures  3  9  through  40B  it  was  necessary 
to  plot  the  profiles  for  each  thermocouple  probe  as  there  were 
too  many  variables  to  produce  a  true  three  dimensional  temper- 
ature profile.   The  x-axis  of  each  graph  is  the  height  of  the 
probe  tip  in  inches  above  the  bed  distributor  plate.   The  y- 
axis  of  each  graph  is  the  temperature  in  °F.   A  probe  grid  is 
provide  for  each  bed  configuration  with  the  probe  number  for 
reference.   In  all  cases  the  temperature  recorded  by  each 
probe  above  the  1.5  inch  height  was  within  the  constant  air 
outlet  temperature  range  specified  above  for  each  data  run. 
All  temperature  data  recorded  can  be  found  in  Appendix  H. 

As  can  be  readily  observed,  the  bed  quickly  attains  a 
constant  temperature  except  for  a  very  small  region  at  the 
bottom  of  the  bed  and  a  narrow  layer  near  the  container  walls . 
The  area  near  the  walls  is  by  assumption  as  the  experimental 
apparatus  was  not  instrumented  to  read  that  area. 

It  appears  that  the  particles  were  heated  by  the  heater 
as  they  transited  the  heater  surface.   Then,  when  they  reached 
the  bottom  of  the  bed,  thermal  equilibrium  with  the  incoming 
cold  air  was  attained  almost  instantaneously.   Except  for  the 
layer  along  the  walls  of  the  bed  and  in  the  small  turbulent 
mixing  region  at  the  bottom  of  the  bed,  thermal  equilibrium 
between  the  particles  and  the  gas  very  nearly  exists  through- 
out the  fluidized  bed. 

Since  the  inside  walls  of  the  bed  container  were  not 
instrumented  with  thermocouples  it  was  not  possible  to 
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Figure  39.   Bed  temperature  vs.  bed  height  (Rc  =  37.5  in  ) 
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Figure  40R.   Bed  temperature  vs.  bed  height  (flc  =  72  in2) 


106 


1  30 


1  30 


120 


Probe  48 


120 


Probe  49 


1  1  o 


1  30 


0.0      0.5      1.0      1.5 


1  1  0 


0.0      0.5      1.0      1.5 


13  0 


120 


Probe  50 


120 


Probe  51 


1  1  o 


1  3  0 


0.0      0.5      1.0 


T 
1.5 


12  0 


1  1  0 


1  30 


Probe  52 


0.0 


1 1 r 

0.5      1.0      1.5 


12  0 


1  1  0 


Probe  54 


1  1  o 


1  3  0 


0.0      0.5      1.0      1.5 


1  2  0 


1  1  0 


1  3  0 


1  20 


1  1  0 


Probe  53 


0.0      0.5      1.0      1.5 


Probe  55 


1 1     i 

0.0      0.5      1.0      1.5 


0.0      0.5      1.0      1.5 

Figure  40B.   Bed  temperature  vs.  bed  height  CRC  =  72  in2) 
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determine  if  there  existed  a  temperature  differential  between 
the  walls  which  did  not  have  heaters  and  the  bed  bulk 
temperature . 

In  order  to  determine  the  effect  of  the  air  velocity  on 
the  temperature  readings  of  the  thermocouple  probes  several 
runs  where  conducted  at  different  air  inlet  velocities  with 
the  probe  tip  placed  six  inches  above  the  top  of  the  distri- 
butor plate.   Small  temperature  fluctuations  of  at  most  1°F 
were  observed  while  the  bed  was  fluidized  and  being  heated 
under  steady  state  conditions.   The  air  flow  was  then  ter- 
minated while  monitoring  the  bed  temperature  readings.   It 
was  observed  that  the  thermocouple  reading  would  initially 
stabilize  at  a  temperature  within  the  range  of  the  previous 
fluctuations  and  did  not  change  for  several  seconds .   A 
gradual  rise  then  would  commence  due  to  the  heat  conduction 
from  the  heaters  in  the  walls.   Since  the  temperature  rise 
did  not  exceed  1°F  at  the  end  of  four  minutes  it  would  seem 
to  indicate  that  the  bed  thermocouple  readings  were  basically 
independent  of  the  fluid  velocity.   It  would  also  tend  to 
support  the  assumption  made  earlier  in  this  thesis  that  the 
particle  temperature  and  gas  temperature  at  any  point  in  the 
bed  can  be  considered  to  be  at  equilibrium. 

The  only  change  in  the  mixing  cell  pattern  within  the 
fluidizing  particle  cell  appeared  to  be  the  production  of 
the  temperature  gradient  at  the  bottom  of  the  bed  and  an 
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increased  bulk  bed  temperature.  This  was  in  comparison  to 
a  steady  state  bulk  bed  temperature  equal  to  the  air  inlet 
temperature  prior  to  energizing  the  plate  heaters. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

The  following  conclusions  may  be  made  from  the  results 
of  this  investigation: 

1.  Upon  fluidization  of  the  bed  an  electrostatic  charge 
will  build  up  within  the  bed.   Although  no  data  correlation 
error  is  used  in  this  thesis,  there  most  probably  is  some 
effect  on  data  results. 

2.  With  a  square  column  bed  the  particle  flow  cell  formed 
is  basically  as  shown  in  Figure  33.   Only  one  major  cell  is 
formed  with  its  axis  of  ascent  roughly  located  at  the  center 
of  the  bed. 

3.  With  an  expanded  rectangular  bed  several  particle 
flow  cells  form  in  a  random  fashion.   There  is  no  pattern 
to  the  overall  cell  coverage  and  the  bed  can  be  termed  a 
bubbling  bed. 

4 .  That  since  the  bulk  bed  temperature  remains  constant 
throughout  the  bed  that  there  must  exist  eddy  cells  within 
the  bed  that  continually  mix  the  bed  particles  as  they  move 
up  through  the  bed.   Otherwise,  a  gradient  would  form  between 
the  heater  walls  and  the  center  of  the  bed. 

5.  That  the  velocity  in  the  center  of  any  cell  was 
greater  than  the  velocity  of  the  particles  being  forced  down 
the  walls.   In  the  expanded  bed  the  particles  moving  down  the 
walls  were  a  mixture  of  the  various  ruptured  cells . 
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6 .  That  mechanically  smooth  f luidization  was  not  taking 
place  at  all  times.  Occasional  instances  of  channeling  were 
present  and  could  be  distinguished  by  small  pertubations 

to  the  bed  surface  that  did  not  change  position  -  a  simmering 
appearance. 

7.  That  the  distributor  plate  used  did  provide  enough 

of  a  bed  pressure  drop  to  allow  fluidization  to  occur  but  that 
the  tendency  for  channeling  or  slug  flow  was  high  if  the  air 
velocity  was  increased  much  past  the  minimum  fluidization 
velocity. 

8.  That  a  temperature  equilibrium  between  the  incoming 
fluidizing  air  and  the  bed  particles  was  reached  almost  in- 
stantaneously in  a  small  mixing  region  of  the  bed  located  at 
the  bottom  of  the  column. 

9.  That  the  only  temperature  gradient  that  existed  in 
the  bed  was  located  at  the  bottom  of  the  bed  where  the  mixing 
region  was  located.   The  remainder  of  the  fluidized  bed  stayed 
at  a  constant  bulk  bed  temperature  irregardless  of  the  number 
of  flow  cells  present. 

10.  That  the  air  velocity  did  not  affect  the  bulk  bed 
temperature  readings.   It  does,  however,  affect  the  depth  of 
the  mixing  region  at  the  bottom  of  the  bed. 

11.  That  energizing  the  plate  heaters  did  not  visibly 
alter  the  particle  cell  flow  pattern  but  that  it  did  produce 
a  gradient  in  the  temperature  profile  of  the  bed  interior. 
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12.  That  the  rectangular  bed  configuration  under  study 
would  trip  into  fluidization  when  the  log  Ap  across  the  bed 
approached  1.25. 

B.   RECOMMENDATIONS  FOR  FURTHER  STUDIES 

Based  on  the  results  of  the  experimentation  to  date, 
research  conducted  in  preparing  this  thesis  and  the  problems 
experienced  in  trying  to  get  the  fluidization  apparatus  to 
work  the  following  areas  of  study  are  recommended  for  future 
thesis  study: 

1.  A  determination  of  what  effect  the  electrostatic 
charge  which  builds  up  during  the  fluidization  process  has 
on  the  heat  transfer  or  flow  cell  formation. 

2 .  A  study  of  the  resulting  changes  in  the  cell  formation 
patterns  when  bed  particles  of  sizes  larger  than  1  mm  are  used 
Particles  in  this  size  range  are  of  much  more  practical  use  in 
the  metallurgical  industry. 

3.  A  study  of  the  flow  cell  patterns  with  beads  smaller 
than  300  microns  which  should  give  enhanced  mechanical 
fluidization . 

4.  A  study  of  the  effect  of  particle  stratification 
(elutriation)  when  particles  of  different  sizes  are  combined 
in  the  bed.   The  combination  can  be  of  varying  sizes/rough- 
ness/sphericity or  any  other  physical  property  of  the  bed 
particles  which  would  promote  the  separation  of  the  particles 
into  distinct  layers  during  fluidization. 
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5.  A  review  of  the  cell  pattern  formation  when  using 
particles  from  various  manufacturers  (same  or  different  bead 
materials)  while  maintaining  a  constant  particle  size. 

6.  Conduct  a  pressure  study  on  the  bed  with  different 
distributor  configurations .   This  would  include  such  designs 
as  porous  plates,  baffles,  calming  chambers,  filter  screens 
or  various  combinations  of  the  aforementioned  designs. 

7.  Study  the  effect  of  changing  the  location  of  the  heat 
source  from  opposite  wall  external  heaters  to  other  configura- 
tional  combinations  such  as  all  four  walls,  adjacent  walls, 
single  wall,  submerged  heater,  etc. 

8.  Conduct  a  velocity  profile  study  similar  to  that  con- 
ducted in  Reference  [43]  .   That  study  used  a  modified  thermis- 
tor anemometer  probe  to  determine  the  flow  pattern  within  a 
circular  column  bed. 

Additionally,  prior  to  conducting  future  research  the 
following  apparatus  modifications  are  suggested: 

1.  Replace  the  left  side  sponge  rubber  gasket  on  the 
movable  wall  with  a  new  gasket. 

2.  Implant  a  pair  of  thermocouples  in  the  unheated  walls 
of  the  fluidized  bed. 

3.  Mount  a  pair  of  thermocouples  on  the  surface  of  each 
strip  heater  to  enable  the  monitoring  of  the  heater  wall 
temperatures . 

4 .  Manufacture  a  stand  for  the  cabling  to  the  thermo- 
couple probes.   The  stand  should  be  of  sufficient  height 
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to  prevent  the  bending  of  a  probe  from  the  weight  of  its 
cable  when  the  probe  is  retracted  from  the  bed. 

5.  Change  the  mounting  method  of  thermocouples  35 
through  38  so  that  they  are  mounted  in  the  backing  plate 
vice  on  the  backing  plate. 

6.  Unless  the  area  of  study  is  to  focus  on  the  electro- 
static charge  built  up  during  f luidization,  provide  some 
means  to  ground  out  this  phenomenon. 

7.  Provide  a  leveling  system  for  the  f luidization  appara- 
tus that  will  ensure  that  the  fluid  flow  direction  out  of 

the  distributor  plate  is  parallel  to  the  bed  container  walls. 

8.  The  probe  taps  on  top  of  the  fluidizing  apparatus 
should  be  bored  to  5/64  inch  for  ease  in  inserting  the  probes 
during  data  runs.   Additional  nylon  bushings  should  also  be 
ordered. 

9.  A  second  flowmeter  of  greater  capacity  than  11.1  scfm 
should  be  installed  in  the  air  supply  line.   This  will  be  re- 
quired to  measure  the  higher  air  velocities  needed  to  fluidize 
the  bed  with  further  expansion. 

10.  A  cooler  air  source  should  be  used  in  future  runs  to 
enhance  the  temperature  differential  at  the  bottom  of  the  bed 
during  f luidization  with  the  wall  heaters  energized. 

11.  A  new  cap  for  the  fill  vent  is  requird  to  prevent  the 
glass  beads  from  getting  between  the  cap  and  plexiglas  thus 
jamming  the  cap. 
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C.   CONCLUDING  REMARKS 

Particle  motion  is  a  complex  function  of  the  particle  and 
gas  properties,  fluidization  velocity,  gas  distributor  and 
geometry  of  the  bed.   While  a  simple  correlation  can  easily 
be  made  it  must  be  tempered  in  its  general  application  by  the 
fact  that  it  will  most  likely  not  give  adequate  allowances 
for  the  way  in  which  the  particles  circulate  and  subsequently 
change  with  changes  in  the  bed  particle  material,  fluidizing 
conditions,  scale  and  design  of  the  fluidization  apparatus. 

The  interpretations  presented  in  this  thesis  are  based 
on  the  premise  of  visible  particle  movement  at  a  given  bed 
height  and  must  be  qualified  by  the  fact  that  they  represent 
a  very  limited  amount  of  data. 

There  is  a  general  opinion  among  researchers  that  the 
bubbles  of  air  passing  through  a  fluidized  bed  grow  in  size 
by  coalescence  of  smaller  bubbles  into  larger  bubbles.   Even- 
tually given  a  deep  enough  bed  the  bubbles  should  coalesce 
into  one  large  bubble  located  in  the  center  of  the  bed  in 
order  to  maintain  the  uniformity  of  the  bed  [43] .   While 
this  opinion  has  been  generated  as  a  result  of  experimenta- 
tion using  circular  beds  it  would  appear  that  it  can  be 
equally  applied  to  rectangular  beds.   No  doubt  had  the 
expanded  rectangular  bed  been  deep  enough,  eventually  it  also 
would  have  conformed  to  the  aforementioned  opinion. 

The  determination  of  particle  flow  patterns  in  a  fluidized 
bed  would  be  of  extreme  value  in  the  establishment  of  a  kinetic 
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or  mass  transfer  model  of  a  fluidized  bed.   The  fact  that 
there  are  a  multitude  of  interrelated  variables  involved, 
unfortunately,  has  made  the  modeling  a  theoretical  tool  vice 
the  needed  operational  scale-up  tool  required  for  industry. 
More  study  will  undoubtedly  be  required  until  a  workable 
industrial  model  can  be  produced. 

It  has  also  been  an  observation  that  while  the  majority 
of  the  fluidization  research  was  conducted  by  American 
researchers  in  the  the  late  forties  and  early  fifties,  that 
the  bulk  of  the  recent  research  has  been  conducted  by  the 
Russians  with  some  notable  exceptions.   This  trend  is  neither 
understood  nor  considered  warranted  by  the  researcher  con- 
sidering the  very  pronounce  heat  transfer  improvement  demon- 
strated through  the  use  of  fluidization. 
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APPENDIX  A:   NONDIMENSIONAL  CORRELATIONS 
FOR  EXTERNAL  WALLS 

These  correlations  are  but  a  sampling  of  the  many 
fluidization  correlations  for  external  walls  that  have 
been  theorized. 

1.   Investigator  -  Baerg  et  al.  [3] 
Correlation  - 


hw  ■  hwmax  -  55  e*P  (  -°-012  {  G  "  °-71  Psb  '  > 


Remarks  -  See  Appendix  C  for  h 


wmax 


2.   Investigator  -  Bartholomew  and  Katz  [5] 
Correlation  - 


w     1.56  +  In  (  Re  Ga"°-S67  -  0.012  ) 


c  ,  G 

Pf 


-0.227  pr°-667  Ga0M2 


Remarks  -  None. 
3.   Investigator  -  Brazelton  [7] 
Correlation  - 


w 


c  .  G' 

Pf 


=  0.72 


D'  G 


-  0.8  7 


(  1  -  £  ) 


-  1.3 


Remarks  -  None. 


117 


4.   Investigator  -  Das  and  Sarkar  [11] 
Correlation  - 


h   D 
w 


=  140  r0'18  Re0'09 


D 


0.2  k 


Remarks  -  None. 
5.   Investigator  -  Dow  and  Jakob  [13] 
Correlation  - 


h   D 
w 


=  0.55 


D 
L 


0.6  5 


D 

d 


0.1  7 


(  1   ~    e    }  ps  cps 


z  pf  cpf 


0.2  5 


0.8  0 


D  G 


Remarks  -  None. 

6.  Investigator  -  Gamson  [17] 
Correlation  - 

st  pr0.6  6  7  =  2>Q  Re  -0.6  9   (  1     _  £   }  -0.3  0 

m 
Remarks  -  None. 

7.  Investigator  -  Lemlich  and  Caldas  [35] 
Correlation  - 

Nu  =  0.055  Re 
Remarks  -  Low  velocity  liquid-f luidization. 
Correlation  - 


St  Pr0-667  =  1.4 


0.7  9 


1 


Ret  e 


Remarks  -  High  velocity  liquid-f luidization 
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8.   Investigator  -  Leva  [36] 
Correlation  - 

Nu  =  0.5 


d   G 

P 


Remarks  -  None. 
9.   Investigator  -  Leva  and  Grummer  [37] 
Correlation  - 


h   =  3000000  k-  D 

w  f   p 


Dp  Gf  n 
V    R 


0.6  0 


Remarks  -  None. 
10.   Investigator  -  Levenspiel  and  Walton  [41] 
Correlation  - 


h   D 

w   r 


0.3 

D       G 

r 

.    Y    J 

.  u    . 

0.5  0 


h   D 
w   r 


=  0.0018 


d      G 

P 

d 
P 

D 

y 

-  1.1  6 


W 


c  ^  G 
Pf 


=  10 


D   G 
r   e 


0»6  5 


h      0.6  (  d   G  ) 
w  P 


0.7  0 


c  c   G 

Pf 
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10.   Remarks  - 


TT  d 


D   = 
r 


6  (  1  -  £  ) 


y  =  (  {  1  +  a2  }U5    -  a3) 


a  =  0.041  (  1  -  e  ) 


D   G 

r 


0.5  o 


11. 


G   =  - 
e 

e 


Investigator  -  Mickley  and  Trilling  [48] 
Correlation  - 


0.0118  p   G 

m 


0.2  6  3 


h   = 
w 


Remarks  -  None. 
12.   Investigator  -  Rao  and  Kaparthi  [56] 
Correlation  - 


Nu  =  0.014  Re1'1 


Remarks  -  Semi-f luidized  bed. 


1  -  £ 


0.i+  0 
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13.   Investigator  -  Richardson  and  Mitson  [57] 
Correlation  - 


Nu  =  119  Pr 


o.i+  o 


c 
ps 

0.2  9 

"k 
s 

pf 

.kf. 

0.0  2 


Re 


N 


D  G 


0.0  7  5 


Remarks  -  Liquid-f luidized  system  with 


N  =  0.02 


3.45  +  — 


14. 


Investigator  -  Rukenshtein  [58] 
Correlation  - 


Nu  =  0.326  Re0-'23  Pr0-33  Ga0'1" 


Remarks  -  None. 
15.   Investigator  -  Toomey  and  Johnstone  [6  2] 
Correlation  - 


log 


h   d 
w   p 


=  0.575  log 


d   u   p  ^ 
P   P   f 


+  0. 130 


Nu  =  3.7  5 


Remarks  -  None. 


d   u  ,-  p^ 
p   mf   f 


log 


u 


u 


mf 


0.1+  7 
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16.   Investigator  -  van  Heerden  et  al.  [66] 
Correlation  - 


Nu  =  0. 58 


B  d   G 

P 


O.k   5 


P     j-     C 

mf   ps 


p,  c  £ 
f   pf 


0.3  6 


0.1  8 


Pr 


0.5  0 


Remarks  -  None. 
17.   Investigator  -  Wen  and  Leva  [72] 
Correlation  - 


Nu  =  0. 16 


c   p   d  «  g0'50 
s   s   p 


0.<+  0 

Re    n 

R 

0.3  6 


Remarks  -  None. 
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APPENDIX  B:   NONDIMENSIONAL  CORRELATIONS 
FOR  IMMERSED  BODIES 

These  correlations  are  but  a  sampling  of  the  many 
fluidization  correlations  for  immersed  bodies  that  have 
been  theorized. 

1.   Investigator  -  Andeen  and  Glickman  [2] 
Correlation  - 


h   Dm 
w   T 


=  900  (  1  -  s  ) 


G  DT  y 


1  0.3  2  6 


Pj.  d 

f   P 


ps  3 


Pr 


0.3 


2. 


Remarks  -  Sometimes  called  the  modified  Vreedenberg 
correlation. 

Investigator  -  Gamson  [17] 

Correlation  - 


St  pr0.6  6  7  =  2>52  Re  -0.8  0 

m 


(  1  -  £  ) 


-  0.3  0 


Remarks  -  Vertical  tubes. 
3.   Investigator  -  Gelperin  and  Einstein  [18] 
Correlation  - 


Nu  =  0.75  Ga0'22 


T 
P 


0.1  it 


Remarks  -  Vertical  tube  bundle. 
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4.   Investigator  -  Gelperin  et  al.  [20] 
Correlation  - 


Nu 


6  (  1  -  £  ) 


=  0.7  3  Re 


0.3  2 


0.3  2 


6  (  1  -  z    ) 


Remarks  -  Horizontal  tube. 
5.   Investigator  -  Glass  and  Harrison  [24] 
Correlation  -  * 


Nu  =  5  7.4 


Dr 


0.9  1 


mf    s 


o.i 


N 


0.3  h 


U 


0.0  8 


P 


6. 


Remarks  -  Staggered  tube  bundle-  N 

Investigator  -  Huntsinger  [29] 
Correlation  - 


u 


mf 


Nu  =  1.4  Reuo'  Pr 


0.5  3   r>~.-  1.6 


Remarks  -  Coil. 
7.   Investigator  -  Kruglikov  [33] 
Correlation  - 

Nu  =  0.6  6  Re0'86  Pr0,29 

Remarks  -  Coil. 


124 


Investigator  -  Miller  and  Logwinuk  [49] 
Correlation  - 


q0.5   v   0.0  7  2   ^   2.1+ 
S  f 


h   =  1.5 


w        d  °'9S  c  -1-6  y°-{ 
P      Pf 

Remarks  -  Vertical  tube  bundle. 
9.   Investigator  -  Petrie  et  al.  [54] 
Correlation  - 


hw  DT 


=  14 


mf 


0.3  3 


Pr 


0.3  3 


Dr 


0.6  e  7 


Remarks  -  Staggered  and  in-line  tube  bundles. 
10.   Investigator  -  Sarkits  et  al.  [60] 
Correlation  - 


0.4  /-.„  0.2  7  r.^.0.3  3 


Nu  =  0.013  3  Reu'*  Ga  u'z  '  Pr 


Remarks  -  Coil. 
11.   Investigator  -  Traber  et  al.  [63] 
Correlation  - 


c 
s 

0.1+5 

D 

0.1  6 

'h 

O 

_cf_ 

d 

P_ 

d 
P 

Nu  =  0.024  Re0'55  Ga0-10  Pr0-33 


Remarks  -  Coil. 


0.1  3 


H 


0.1  6 
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12.   Investigator  -  Vreedenberg  [69] 
Correlation  - 


Nu  =  a 


G  c   V    ~ 

mf   mf 


Remarks  -  Tube  position  varies. 
13.   Investigator  -  Vreedenberg  [70] 
Correlation  - 


0.3  5 


Nu  =  0.66 


G  D  p   (  1  -  £  ) 
s 


pf  y  e 


0.1+  i+ 


Pr 


0.3  0 


Remarks  -  Horizontal  tube  with- 


Correlation  - 


d   G  o 

P    s 


y  P. 


<  2050 


Nu  =  420 


G  D  y 


pf  V  ps  9 


Remarks  -  Horizontal  tube  with- 


0.3  0 


Pr 


0.3  0 


d   G  p 
p     s 


y  p. 


>  2550 


Remarks  -  When 


d   G  p 
P    s 


y  p. 


is  greater  than  2050  but 


less  than  2550  the  value  of  Nu  will  be 
the  average  of  the  two  Nu  correlations 
above . 
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14.   Investigator  -  Wender  and  Cooper  [74] 
Correlation  - 


Nu  =  0.0119  (  1  -  z    ) 


f          pf 

O.k  3 

C 

ps 

pf. 

0.8 

V 

.pf. 

k 
s 

0.6  6 


Re0-23  C 


Remarks  -  Vertical  tube. 
15.   Investigator  -  Zenz  and  Othmer  [80] 
Correlation  - 


Nu  =  14 


-I  0.3  3 


mf 


Pr 


0.3  3 


Dr 


0.6  7 


Remarks  -  Horizontal  tube. 


R 
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APPENDIX  C:   NONDIMENSIONAL  CORRELATIONS  FOR  THE 
MAXIMUM  HEAT  TRANSFER  COEFFICIENT 

These  correlations  are  but  a  sampling  of  the  many 
fluidization  correlations  for  maximum  heat  transfer  coef- 
ficients that  have  been  theorized. 

1.   Investigator  -  Baerg  et  al.  [3] 
Correlation  - 


h      =49  log 
wmax        ^ 


0.00037  p  , 
sb 


Remarks  -  None. 
2.   Investigator  -  Chechetkin  [8] 
Correlation  - 


Nu 


max 


0.0017  Re    °'8  Pr0-1* 
opt 

d  °'69 
P 


Remarks  -  Re     =  0.209  Ga0,52. 
opt 

3.   Investigator  -  Chekansky  et  al.  [9] 
Correlation  - 


h      =  28.2  p 
wmax         s 


0.2 


k   °'6   d 

f       P 


-  0.3  6 


A 

d 
P_ 

0.0  it 

'DT    ' 

_D20. 

-  0.1  2 


Remarks  -  Staggered  tube  bundle 
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4.   Investigator  -  Gelperin  et  al.  [19] 
Correlation  - 


Nu     =0.74  Ga°*22 
max 


Dr 


D. 


T  /         T 
1    -   1  +  

y  V    z  +  Dr 


0.2  5 


Remarks  -  Staggered  horizontal  tube  bundle. 
5.   Investigator  -  Gelperin  et  al.  [21] 
Correlation  - 


Nu     =  0.6  4  Ga 
max 


0.2  2 


P 


0.0  9 


Remarks  -  Vertical  tube  bundle. 
6.   Investigator  -  Gelperin  et  al.  [22] 
Correlation  - 


Nu     =0.79  Ga0,22 
max 


T 
1 

y 


0.2  5 


Remarks  -  In-line  horizontal  tube  bundle. 
7.   Investigator  -  Jacob  and  Osberg  [30] 
Correlation  - 


h     =  h   (  1  -  e  )  (  1  -  exp  {  -pkp  }  ) 
wmax     o  r     r  f 


Remarks  -  h   and  p  determined  by  experimentation 
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8.   Investigator  -  Maskayev  and  Baskakov  [44] 
Correlation  - 


Nu     =  0.21  Ga0-32 
max 


9. 


Remarks  -  Vertical  cylinder  in  coarse  particle  bed 
with  1.4xl05  <  Ga  <  3xl08. 

Investigator  -  Pillai  [55] 

Correlation  - 


Nu     =  0.365  e0'82  Ga0'22 
max 


Remarks  -  Sphere  with  Bi  <  0.04 
10.   Investigator  -  Ruckenstein  [58] 
Correlation  - 


Nu     =  Re   .  °-k2"    Ga0'1"    PrM3 


max 


'opt 


Remarks  -  Re     =0.09  Ga0,58 
opt 

11.   Investigator  -  Sarkits  [60] 
Correlation  - 


Nu     =  0.0087  Ga0-1*2  Pr0'33 
max 


c 

ps 

QMS 

D 

0.1  6 

"h 

o 

pf 

U                 J 

d 
P_ 

d 

P_ 

0A  5 


Remarks  -  Re     =  0.2  Ga°'D  (laminar  region). 

opt  3 

Correlation  - 


Nu     =0.019  Ga0'5  Pr0,3  3 
max 


c 

ps 

0.1 

D 

0.1  3 

rH     1 
o 

c     xz 
Pf 

— 

d 
P 

d 

P 

0.1  6 


Remarks  -  Re     =  0.66  Ga0-5  (turbulent  region). 

opt  3 


130 


12.   Investigator  -  Traber  et  al.  [63] 
Correlation  - 


Nu     =  0.021  Ga0'1*  Pr0-33 
max 


D 

0.1  3 

"H    1 

O 

d 
P 

d 
P 

L        _J 

-      —j 

0.1  6 


Remarks  -  Re   t  =  0.55  Ga0,5. 
opt 

13.   Investigator  -  Varygin  and  Martyushin  [67] 
Correlation  - 


Nu     =  0.86  Ga 
max 


0.2 


Remarks  -  Re   .  =  0.118  Ga0,5. 
opt 

14.   Investigator  -  Zabrodsky  [78] 
Correlation  - 


h      =  33.7  o  °-2  k  °'G  d  -0-36 
wmax         s     f     p 


Remarks  -  None. 
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APPENDIX  D:   GENERAL  EQUATIONS  AND  UNITS 

This  appendix  will  be  broken  down  into  two  sections. 
The  first  section  will  cover  the  general  nondimensional  heat 
transfer/fluids  numbers  used  in  this  thesis.   The  second 
section  will  annotate  the  units  to  be  used  with  each  equation 
that  is  numbered  within  the  main  text.   This  requirement  was 
necessitated  by  the  use  of  different  units  (same  variable 
nomenclature)  among  the  various  researchers  cited  in  the 
thesis  text. 
A.   GENERAL  NONDIMENSIONAL  HEAT  TRANSFER/FLUIDS  NUMBERS 

1.   Reynolds  Number  (Re) 


d   u   p_ 
p  o  wf 


Re  = 


2.   Galileo  Number  (Ga) 


V  pf  (  ps  "  pf  }  g 


Ga  = 


u2 


3.   Prandtl  Number  (Pr) 


c  -  u 
pf  M 


Pr  =  

k 


f 
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4.   Nusselt  Number  (Nu) 


h   d 
w  p 

Nu  =  

kf 


5.   Stanton  Number  (St) 


h 

w 


St  = 


cpf  pf  Uo 


6.   Biot  Number  (Bi) 


h   r 
w   sp 

Bi  =  

k 
sp 


B.   EQUATION  UNITS 

1.   Equation  (3.1) 


4pe 

- 

lb/ft2 

L 
e 

- 

ft 

e 

- 

dimensionless 

ps 

- 

lbm/ft3 

pf 

- 

lbm/ft3 

2.   Equation  (3.2) 

mf 

- 

dimensionless 

V 
P 

- 

ft3 

Yf 

- 

lb/ft3 

Lmf 

- 

ft 

A 
c 

- 

ft2 

ps 

- 

lbm/ft3 
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pf   -  lbm/ft 


3.   Equation  (3.3) 


G  ,   - 
mf 

lbm/(hr-ft2) 

D 

P 

in 

y 

centipoise 

ps   " 

lbm/ft3 

pf   " 

lbm/ft3 

Equation  (3.4) 

5 

P 

in 

D. 

l 

in 

x. 

l 

dimensionless 

Equation  (3.5) 

Remf  " 

dimensionless 

Gmf   " 

lbm/(hr-ft2) 

5 

p 

in 

y 

centipoise 

Equation  (3.6) 

Gmf   -  lbm/(hr-ft2) 

Equation  (3.7) 

n    -  dimensionless 

Gf    -  lbm/(hr-ft2) 

G    -  lbm/(hr-ft2) 

Equation  (3.8) 

r     -  dimensionless 

m"    -  dimensionless 

Gf    -  lbm/(hr-ft2) 
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G  -   -  lbm/(hr-ft2) 
mf 

9.   Equation  (3.9) 


q     -  BTU/hr 

h  ._  -  BTU/(hr-ft2-°F) 
wiL 

A     -  in2 
m 

T     -  °F 
w 

TK    -  °F 
b 


10.   Equation  (3.10) 


2   0 


h  ._  -  BTU/(hr-ft--°F) 
wiL 


k     -  BTU/(hr-ft-°F) 
m 


p     -  lbm/ft 
m 


3 


c     -  BTU/(lbm-°F) 
m 

t    -  sec 

11.   Equation  (3.11) 

h     -  BTU/(hr-ft2-°F) 
w 

k     -  BTU/(hr-ft-°F) 
m 


P    -  lbm/ft 

m 


3 


c  -  BTU/(lbm-°F) 
m 

u  -  ft/hr 
s 

LH  -  ft 

12.   Equation  (3. 12) 

h  T  -  BTU/(hr-ft2-°F) 

k  -  BTU/(hr-ft-°F) 
m 

cm  -  BTU/ ( lbm- ° F ) 

s  -  1/sec 
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APPENDIX  E:   EQUIPMENT  LISTING 


Device 

Manufacturer 

Model  # 

Serial  # 

Temperature 
Recorder 

Newport 

267B-TF2 

8500392 

Temperature 
Recorder 

Newport 

267B-TF2 

8500393 

Thermocouple 

Omega 
Engineering 

ANSI  T 
Type  S 
Sheathed 



Water 
Manometer 
(30  in) 

Meriam 
Instrument 

Type  W 
33KA35 

A69513 

Water 
Manometer 
(60  in) 

Meriam 
Instrument 

Type  W 
A324 

44F554 

Flowmeter 

Fischer  and 
Porter 

10A3565A 

6908A2165A1 

Variac 

General  Radio 

Type  V10 



Air  Blower 

Clements 
National 

G12 

114986 

Regulated 
Power  Supply 

Lambda 
Electronics 

LK  345A  FM 

C  1821 

Thermocouple 
Calibration 

Rosemount 
Engineering 

Bridge-920A 
Bath-913A 

110 

Strip 
Heaters 

Wat low 
Electric 

250  Watts 

8044 

Weight  Scale 

Mettler 
Instrument 

Type  H15 

119923 

Optical 
Microscope 

Bausch  &  Lomb 
Instrument 

31-32-45 

B-3539  RS 
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APPENDIX  F:   PARTICLE  ANALYSIS 

The  particle  specifications  provided  by  the  manufacturer, 
Potters  Industries,  are  as  follows: 

1.  Physical  Properties 

a.  Size  -  Minimum  90%  by  weight  in  range  on  smaller 
U.S.  Sieve,  maximum  of  10%  by  weight  on  larger  U.S.  Sieve. 
Tested  in  accordance  with  ASTM-D-1214-58  and  MIL  Spec. 
G-9954A. 

b.  Free  Iron  Content  -  Maximum  of  0.1%  by  weight. 

c.  Chemical  Content  -  Soda-lime  silica  glass.   The 
solid  glass  spheres  contain  no  free  silica. 

d.  Specific  Gravity  -  2.45-2.55. 

e.  Hardness  -  DPH  lOOg  load-515Kg/mm2 . 

f.  Color  -  Clear  and  colorless,  free  from  surface 
films . 

g.  Roundness  -  Minimum  percent  round  by  roundometer 
is  90%.   Determined  in  accordance  with  ASTM-D-1155-53  or 
MIL  Spec.  G-9954A. 

h.   Broken  or  Angular  Particles  -  Maximum  of  3%  by 
count.   Tested  in  accordance  with  MIL  Spec.  G-9954A. 

i.   Air  Inclusions  -  Not  more  than  10%  of  the  beads 
shall  show  air  inclusions  of  more  than  25%  of  their  surface. 
Tested  microscopically  immersed  in  1.5  refractive  index 
fluid. 

j.   Average  Settled  Bed  Density  -  96  lbm/ft3. 

k.   Average  Particle  Density  -  156  lbm/ft3. 

1.   Coefficient  of  Friction  (Static)  -  0.9-1.0. 

2.  Thermal  Properties  @  70°F 

a.   Specific  Heat  -  0.180  BTU/ (Ibm- °F) . 
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b.   Thermal  Conductivity  -  0.515  BTU/ (hr-f t- °F) . 

3.  Electrical  Properties  @  77°F 

a.  Dielectric  Constant  (1  Kc)  -  7.6. 

b.  D.C.  Volume  Resistivity  -  6. 5X10 12  Ohm-cm. 

4 .  Size  Analysis  (Potters  P-020  Designation) 

a.  Basic  Sieve  Analysis  -  U.S.  Sieve  35  to  50. 

b.  Average  Particle  Diameter  -5   =  0.0122  inch 
(310  microns) .  p 

c.  Maximum  Percentage  True  Spheres  -  90%. 
Photographs  of  two  groups  of  the  particles  are  shown  in 

Figures  41  and  42.   Upon  microscopic  examination  the  beads 
appeared  extremely  spherical  in  nature  with  some  air 
inclusions  present  in  all  beads.   The  inclusions  appear  as 
black  areas  on  the  photographs.   There  was  little  variation 
in  particle  diameters  between  the  two  sample  groups.   The 
diameter  for  each  group  was  taken  as  the  arithmetic  mean  of 
microscopic  measurements  on  a  representative  sample  of  50-60 
beads  in  each  group.   Group  I  had  an  average  particle 
diameter  of  0.0111  in  while  Group  II  had  an  average  particle 
diameter  of  0.0110  in.   For  the  purpose  of  calculations  in 
this  thesis  an  average  particle  diameter  (D  )  of  0.011  in 
was  used. 

A  basic  sieve  analysis  was  performed  on  two  samples  of 
known  volume  and  weight  to  verify  the  manufacturer's 
specifications.   The  sieves  where  U.S.  Tyler  Company  stand- 
ard screen  mesh  of  48  and  60  mesh  size.   Less  than  five 
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Figure  41.   Mi crophotogr aph  of  beads-Group  I 


Figure  42.   Mi crophotogr aph  of  beads-Group  II. 
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percent  by  weight  of  the  sample  particles  was  larger  than 
the  48  mesh  screen  and  less  than  one  percent  by  weight  of 
the  sample  was  smaller  than  60  mesh.   The  Tyler  screen  scale 
of  48  (U.S.  Sieve  No.  50)  has  a  sieve  opening  of  approximate- 
ly 0.0117  in  and  the  Tyler  screen  scale  of  60  (U.S.  Sieve 
No.  60)  has  a  sieve  opening  of  approximately  0.0098  in. 

Water  displacement  of  a  known  weight  of  beads  was  used 
to  determine  the  average  particle  density.   The  value 
obtained  was  154.6  lbm/ft3  which  is  considered  to  be  within 
acceptable  agreement  of  the  manufacturer's  specifications. 
Likewise,  the  average  settled  bed  density  was  obtained  by 
measuring  the  volume  occupied  by  a  packed  bed  of  beads  of 
known  weight.   The  value  obtained,  93.2  lbm/ft3,  is  also 
within  acceptable  limits. 
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APPENDIX  G:   EXPERIMENTAL  UNCERTAINTY  ANALYSIS 

The  source  of  error  in  the  results  presented  is  consid- 
ered to  be  due  to  instrument  precision  and  accuracy, 
inaccuracies  in  geometrical  measurements  and  inaccuracies 


in  data  recording.   The  following  is  a  breakdown  of  the 
experimental  uncertainties. 

1.   Distance  Measurements 

The  tape  attached  to  the  side  of  the  fluidizing 
apparatus  has  a  minimum  scale  division  of  0.25  in.   Distance 
readings  were  made  with  this  scale  to  the  nearest  0.0625  in 
for  both  the  bed  height  and  width.   While  it  was  possible  to 
get  static  readings  with  a  precision  of  ±0.5%,  the  expanded 
bed  heights  were  more  on  the  order  of  ±2%  with  a  worst  case 
precision  of  ±4%. 

The  depth  measurement  of  the  thermocouple  probes  was 
based  on  measuring  the  removal  distance  of  the  probe  from 
its  tap  at  the  top  of  the  apparatus.   The  smallest  scale 
division  for  this  withdrawal  measurement  was  0.0625  in. 
However,  since  the  probes  were  not  rigid  structures,  the 
bending  in  the  probe  lengths  resulted  in  precision  estimates 
of  ±2%. 

12  .   Particle  Diameter  Measurements 
The  optical  microscope  used  to  measure  bed  particle 
diameters  in  the  two  samplings  of  beads  has  a  minimum  scale 
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division  of  0.0001  in  which  allowed  reliable  measurements  up 
to  0.00005  in  at  100X  magnification.   The  diameter  precision 
based  on  scale  readings  is  ±0.5%.   However,  operator  error 
in  reading  the  diameters  places  the  precision  estimate  more 
realistically  at  ±1%.   The  standard  deviation  in  particle 
diameter  for  the  two  bead  groups  was  0.00042  in  and  0.00045 
in. 

3 .  Particle  Density  Measurements 

The  scale  used  to  measure  the  particle  bead  weights 
has  a  minimum  division  of  1  mg  which  equates  to  a  precision 
of  ±0.6%.   The  minimum  volume  scale  was  1  ml  which  equates 
to  a  precision  of  ±1%.   Due  to  spillage  in  filling  the 
measuring  containers  with  water  and  beads  in  addition  to  the 
losses  during  the  experimental  runs  it  is  estimated  that  the 
error  in  density  measurements  is  closer  to  x3%  overall. 

4 .  Temperature  Measurements 

Although  the  temperature  recorder  gives  readouts  in 
tenths  of  degrees,  after  taking  into  consideration  the 
fluctuations  experienced  during  fluidization  and  round  off 
errors  the  overall  error  in  temperature  measurements  is 
considered  to  be  ±1°  F. 

5.  Pressure  Measurements 

The  minimum  manometer  scale  division  was  0.1  in  H2O 
with  readings  of  0.05  in  H20  possible  for  a  full  scale 
precision  of  ±0.2%.   While  this  precision  was  possible  during 
static  bed  conditions  it  was  not  possible  to  be  as  accurate 
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once  fluidization  occurred  as  the  pressure  level  would 
fluctuate  within  the  manometer  column.   For  the  worst  case 
of  the  pressure  measurement  in  the  column,  the  pressure 
fluctuated  ±0.3  in  H2O  over  the  maximum  pressure  recorded 
for  a  recorded  error  accuracy  of  ±2%. 

6 .  Flow  Measurements 

The  minimum  flowmeter  division  was  1%  of  11.1  scfm 
or  0.111  scfm.   This  gives  a  scale  precision  of  ±1.2%  for 
the  range  measured.   Precision  in  reading  the  float  position 
especially  with  fluctuations  during  fluidization  caused 
additional  inaccuracies.   As  a  result,  the  overall  flow 
measurement  was  probably  within  ±5%. 

7.  Miscellaneous  Measurements/Calculations 

The  following  uncertainties  were  determined  for 
various  calculations/measurements  required  for  data  reduction 
by  the  method  of  Kline  and  McClintock. 

a.  A   -  ±1.5%. 

c 

b.  Gf  -  ±5. 7%. 


c.  pf  -  ±2.4%. 

d .  Re  -  ±5.9%. 

e.  y   -  ±1.0%. 


143 


APPENDIX  H:   EXPERIMENTAL  DATA 

The  following  temperature  and  pressure  data  was  recorded 
during  experimental  runs  of  the  fluidized  bed  apparatus. 

1.  Pressure  Data  -  Run  #1 

a.  Static  Bed  Height  -  12.5  inches. 

b.  Static  Bed  Width  -  6.25  inches. 

c.  Air  Inlet  Cross-Sectional  Area  -  37.5  in2. 

d.  Bed  Pressure/Gas  Mass  Flow  Data  -  See  Table  1. 

e.  Comments  - 

(1)  Channeling  observed  to  start  at  the  edges 
of  the  bed  at  an  air  flow  of  8.88  scfm. 

(2)  Fluidization  of  the  bed  commenced  at  an  air 
flow  of  9.66  scfm. 

(3)  Bed  expansion  during  fluidization  was 
0.125-0.375  inch. 

2.  Pressure  Data  -  Run  #2 

a.  Static  Bed  Height  -  12.5  inches. 

b.  Static  Bed  Width  -  8.0  inches. 

c.  Air  Inlet  Cross-Sectional  Area  -  48.0  in2. 

d.  Bed  Pressure/Gas  Mass  Flow  Data  -  See  Table  2. 

e.  Comments  - 

(1)  Channeling  observed  to  start  at  the  edges 
of  the  bed  at  an  air  flow  of  9.44  scfm. 

(2)  Fluidization  of  the  bed  commenced  at  an  air 
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flow  of  11.32  scfm. 

(3)   Bed  expansion  during  fluidization  was 
0. 125-0. 375  inch. 

3.  Pressure  Data  -  Run  #3 

a.  Static  Bed  Height  -  12.5  inches. 

b.  Static  Bed  Width  -  10.0  inches. 

c.  Air  Inlet  Cross-Sectional  Area  -  60.0  in2. 

d.  Bed  Pressure/Gas  Mass  Flow  Data  -  See  Table  3. 

e.  Comments  - 

(1)  Channeling  observed  at  the  corners  of  the 
bed  at  an  air  flow  of  10.55  scfm. 

(2)  Channeling  observed  along  the  walls  and 
occasionally  in  the  bed  at  an  air  flow  of  11.66  scfm. 

(3)  Fluidization  of  the  bed  commenced  at  an  air 
flow  of  12.77  scfm. 

(4)  Bed  expansion  during  fluidization  was  0.25 
inch. 

4.  Pressure  Data  -  Run  #4 

a.  Static  Bed  Height  -  12.5  inches. 

b.  Static  Bed  Width  -  12.0  inches. 

c.  Air  Inlet  Cross-Sectional  Area  -  72.0  in2. 

d.  Bed  Pressure/Gas  Mass  Flow  Data  -  See  Table  4. 

e.  Comments  - 

(1)  Channeling  observed  at  the  corners  of  the 
bed  at  an  air  flow  of  11.66  scfm. 

(2)  Channeling  commenced  in  the  main  part  of 


145 


the  bed  at  an  air  flow  of  12.77  scfm. 

(3)  Fluidization  of  the  bed  commenced  at  an  air 
flow  of  15. 00  scfm. 

(4)  Bed  expansion  during  fluidization  was  0.25 
inch. 

5 .  Temperature  Data  -  Run  #1 

a.  Static  Bed  Height  -  12.5  inches. 

b.  Static  Bed  Width  -  6.25  inches. 

c.  Air  Inlet  Cross-Sectional  Area  -  37.5  in2. 

d.  Average  Air  Inlet  Temperature  -  108°  F  ±2.5°  F. 

e.  Constant  Air  Outlet  Air  Temperature  -  125°  F 
±0.4°  F. 

f.  Temperature  Distribution  Data  -  See  Table  5. 

6 .  Temperature  Data  -  Run  #2 

a.  Static  Bed  Height  -  12.5  inches. 

b.  Static  Bed  Width  -  12.0  inches. 

c.  Air  Inlet  Cross-Sectional  Area  -  72.0  in2. 

d.  Average  Air  Inlet  Temperature  -  116°  F  ±2.0°  F. 

e.  Constant  Air  Outlet  Air  Temperature  -  129.5°  F 
±0. 7°  F. 

f.  Temperature  Distribution  Data  -  See  Table  6. 
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